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ABSTRACT 

Photoelastici ty gives the inf oriBation of ia - ty ) in the form of 

1 2 

fringes. Fringe ordering is one of the most difficult steps in 
processing the fringe patterns from any optical technique. Automated 
ordering of the fringes remains as one of the most challenging tasks 
even today. In this thesis, a new methodology is evolved for automatic 
fringe ordering using the technique of D.I.P. The new methodology has 
three steps <1) fringe skeletonisation <2) fringe tracing and C3) fringe 
gradient identification. Fringe skeleton is obtained using the algorithm 
of Ramesh and Pramod. Fringe tracing algorithm is newly developed and it 
identifies the pixels forming a continuous fringe. The algorithm uses a 
3X3 mask and a tree structure for tracing out the fringes. Intensity 
variation over the fringe field is used to identify the fringe gradient 
direction. It is established in this thesis that the intensity variation 
over the fringe field is not trivial .Complete automation of the 
procedure has been implemented on a PC based system for a class of 
fringe patterns such as those occurring in circular disc under diametral 
compression, A semi-automatic procedure has been implemented for the 
fringe field occurring in the neighborhood of a crack. 

Most of the practical structures of engineering importance have 
finite dimensions and the cracks in them may be located in zones of 

stress concentration and the cracks may be large enough so that the 

crack— tip may be closer to a boundary. Multi parameter stress field 
equations are needed to evaluate the stress field. Among the 
mul ti-parameter stress field equations, the 'r.:_,^ci'cions proposed by Atiuri 
and Kobayashi are the simplest. Several typographical errors in the 
reported equations are corrected in this thesis. Using these equations 
and following the methodology of Sanford, an overdeterministic least 
square method is developed to extract the stress field parameters from 
the isochroma tic fringes. The application of this methodology to a 
practical problem of the mixed mode analysis of a crack at the tensile 

fillet of a spur gear is discussed- Results show that proper choice of 

the number of parameters can model the stress field in the vicinity, as 
well as the far field, accurately. 
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CHAPTER 1 

INTRODUCTION 


1.1 Introduction 

When a ray of polarised light is incident on the surface of a 
model made of birefringent material, it splits up into two 
components, one of which obeys the Snell's law the other < 
extraordinary component ) violates the Snell's law and moves at a 
slower speed. This leads to a phase lag between the two rays as they 
emerge out of the model, creating a retardation pattern. 

It has been found out experimentally that this phase lag is 
directly proportional to the difference of the principal stress and 
the stress-optic law relates this as, 

N F 

— = i O' — o ) < 1 . 1 ) 

t * ^ 

where <N/mm/f ringe) is the material fringe constant, N is the 
fringe order and t is the model thickness. The contours 
corresponding to are known as isochromatics . 

If a plane polarised light is used, then one observes two sets 
of contours namely isochromatics (contours of maximum shear stress) 
and isoclinics (contours of principal stress orientation). If a 
circularly polarised light is used, then one observes only 
isochromatics . 

An experiment in photoelastici ty involves the following steps : 

1. Image acquisition. 

2. Fringe extraction and ordering 

3. Data collection and its analysis. 

The collection of data involves the determination of the fringe 
order and the positional coordinates of various data points from the 
fringe field. 
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1 . 2 Pr eseni wor k 

One of the main advantages of a technique like Photoelasticity 
is that one gets the information related to the stress field rather 
than a point information as obtained in strain gauge. The 
instrumentation in strain gauge has developed to an extent that one 
can directly read the values from the meter as strains. On the other 
hand, one obtains the result in Photoelasticity in the form of 
fringes. Though the fringes must appear as lines theoretically, they 
appear as broad bands in view of the limited resolution of the 
recording medium. The breadth of the band is related to the local 
stress gradient. 

In order to obtain quantitative information, one has to 
identify first the exact fringe location and then find the fringe 
order. With the development of Digital Image Processing, now 
sophisticated algorithms are available to extract the fringe 
skeletons precisely. 

Though the fringe skeleton can be extracted, there is no 
algorithm available for automatic fringe ordering. Fringe ordering 
is one of the most difficult steps in processing the fringe patterns 
from any optical technique. One needs prior knowledge to do this and 
automated ordering of the fringes remains as one of the most 
challenging tasks even today. 

In this thesis, a new methodology is evolved for automating the 
fringe ordering in photoelastic experiment. Using the fringe 
thinning algorithm by Ramesh and Pramod ClI the fringe skeleton is 
extracted. An algorithm is developed and implemented for fringe 
recognition and tracing. Complete automation of the procedure is 
achieved for certain class of problems. 

The present work is a step towards complete automation for a 
Photoelastic experiment. A comprehensive software can be developed 
as an extension which will be able to interpret thf' fringe patterns 
and give the results directly in the form of ^tresses for the 
photoelastic experiment. 
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This will make the photoeiastic technique accessible to a much 
broader spectrum of users. Combined with the inherent advantages of 
photoelasticity, i.e. the availability of the whole field 
information, this technique would therefore lead to wider 
application of photoelasticity in the industry. 

The use of Fracture Mechanics principles to design structures 
is gaining importance in the recent years and one of the important 
parameter to be evaluated is the Stress Intensity Factor (S.I.F). 
The S.I.F can be evaluated by analysing the isochromatic fringe 
field in the neighborhood of the crack-tip. One must use appropriate 
equations, which closely model the stress field for data reduction. 

Most of the practical structures of engineering importance have 
finite dimensions and the cracks in them may be located in zones of 
stress concentration and the cracks may be large enough so that the 
crack-tip may be closer to a boundary. The basic stress field 
equations of fracture mechanics are developed for cracks lying in a 
body of infinite dimensions( i .e« , the crack-tips are far away from 
the boundaries) subjected to an uniform stress field. When these 
equations are attempted to be applied for solving real life 
problems, one gets erroneous results. The recognition of this fact 
has lead the researchers to account for this effect and 
multi-parameter stress field equations have been developed in the 
past two decades. 

Among the various stress field equations, the use of equations 
proposed by Atluri and Kobayashi C23 is selected for their 
compactness and ease of programming. Though the equations proposed 
by Atluri and Kobayashi are recent ones and elegant to use, 
unfortunately there were several typographical errors in the 
reported equation. The errors were identified when the reconstructed 
fringe field did not have symmetry in the case of Mode II. These 
errors are corrected in the present thesis. Following the lines of 
Sanford £33 an overdeterministic scheme is evolved to solve for the 
mixed mode parameters. The program is validated by solving a Mode-I 
crack problem. Use of the above methodology for solving the problem 
of finding the stress field in the neighborhood of a crack in the 
tensile fillet of a spur gear is discussed towards the end. 
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1.3 T^iesls layout 

Chapter 2 deals with the new methodology for automatically 
tracing the fringes. The fringe thinning algorithm scans the image 
along particular directions and marks all the fringe points, however 
it cannot recognize a particular fringe as a continuous one since, 
in the computer representation, the fringes are just a group of 
points without any particular ordering. The first step towards 
automatic fringe ordering lies in the computer recognition of the 
fringes. Starting with a simple 3X3 mask, its performance is 
evaluated for various test images and also for an actual fringe 
pattern. Improvements are made in this and finally a complex fringe 
tracing algorithm is presented which improves upon the shortcomings 
of the previous algorithms. 

Various aspects involved in the fringe ordering procedure in 
photoelasticity are discussed in Chapter 3. The importance of the 
intensity variation in the fringe field to identify fringe gradient 
direction in photoelasticity is established. An algorithm is 
presented for complete automation of the fringe ordering procedure 
for a circular disk under diametral compression and an interactive 
approach is proposed for automation of the fringe ordering for the 
mixed mode crack problem. 

Chapter 4 deals with the least square procedure for the 
analysis of the data collected to evaluate the crack-tip stress 
field parameters from the photoelastic experiment. The need for 
considering the higher order terms to model the stress field in the 
neighborhood of long cracks is established. Using the Atluri's 
series solution C21, an example problem of Mode I crack is analysed 
and the fringe field is reconstructed using computer graphics for 
comparison with the experimental patterns. The typographical errors 
of the reported equations in the original referenceC2] are also 
brought out. 

In Chapter 5 application of the above procedure to mixed mode 
problem is discussed. The problem of evaluating the stress field 
parameters in the neighborhood of a crack at the tensile fillet of a 
spur gear is studied. 
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One of the major tasks completed in the course of the present 
investigation is the conversion of various software modules 
developed around Professional Image Processing < PIP ) system over 
the last four years to the newly acquired MVP~AT system. Software 
modules for image acquisition, image filtering and cleaning , fringe 
skeletonisation using the algorithm of Ramesh and Pramod C13 have 
been implemented. Appendix section lists out the prominent 
modifications required for this. 
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CHAPTER E 

AUTOMATED FRINGE RECOGNITION BY FRINGE TRACING 


2.11 rjiroductl on 

The photoelastic -fringes generally appear in the -form of bands 
and not as thin lines- Identification of actual fringe skeleton is 
an important task. The accuracy in the data collection is subject to 
human judgment in the conventional analysis. 

However, with the development of Digital Image Processing, now 
sophisticated algorithms are available to extract the fringe 
skeletons precisely. The fringe thinning algorithm scans the image 
along particular directions and marks all the fringe points, however 
it cannot recognize a particular fringe as a continuous one since, 
in the computer representation, the fringes are just a group of 
points without any particular ordering. The first step towards 
automatic fringe ordering lies in the fringe recognition algorithm. 
In this chapter, a fringe tracing algorithm is proposed and its 
performance is accessed for various test images and actual 
photoelastic fringe patterns. 

2.2 System Conf igur atlon 

The image processing system used consists of a CCD camera (TH 
560 PULNIX) that has a pixel resolution of 512 X 512 pixels and 
digitizes the image at video rates. 

The camera is connected to a PC— 386 based colour 
image-processing system equipped with an image processing card 
<HVP-AT, Matrox corporation) connected to a high resolution video 
monitor as shown in Fig. 2.1 The system has in addition to PC memory, 
four frame buffers, each of 512 X 512 pixels X 8 bits, to store the 
image being processed and the intermediate results. 
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2.3 Fringe Skeleton Extraction 

In Digital Image Proceesing (D.I.P), the image is identi-fied as 
an assembly of picture elements abbreviated as pixels. Depending on 
the intensity of light transmitted/reflected each pixel is assigned 
a number between 0 (pitch dark) and 255 (pure white). Now the image 
is available as a matrix of numbers which can be handled by a 
digital computer for further processing. Due to the nature of 
formation of fringes, the intensity of pixels in the fringe band 
does not remain constant but there exists a variation. Though the 
variation is not resolvable to a fine degree by a human eye, the 
present day CCD cameras can easily recognise it and quantise it for 
further processing. 

The various DIP algorithms reported in the literature can be 
classified into two categoriesC41 . In the first category the fringe 
field is identified as a binary image and the fringe skeletons are 
obtained using algorithms that were primarily developed for optical 
character recognition. The algorithms of MQller and Saackel C53, 
Seguchi et.al. C61 and the recent one of Chen and Taylor C71 come 
under this category. In the other category, intensity variation 
within a fringe is used in one way or other in devising algorithms 
for fringe skeleton determination. The algorithms of Yatagai et.al. 
C43, Gillies C93, Umezaki et.al. C103, and Ramesh and Pramod C13 
come under this category. 

Utilisation of intensity variation, within a fringe, to detect 
fringe skeleton was first reported by Yatagai et.al. C83. They were 
essentially concerned with the detection of fringe maxima (in most 
systems 0 is pitch black and 255 is pure white) to identify fringe 
skeletons from the bright bands. Umezaki et.al. C103 reported another 
algorithm for fringe skeleton detection using fringe minima as a 
criterion and concentrated in extracting fringe skeletons from dark 
bands. 
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£.3.1 Algor Itlufi of Rainesh and Prai»od II) 

Ramesh et.al. till showed for the first time that fringe edge 
detection followed by fringe skeleton identification using a 
minimum intensity criterion greatly minimises the noise reported in 
other algorithms. 

However, the algorithm is limited to processing either 
horizontal or vertical fringes. Ramesh and Pramod £11 improved the 
earlier algorithm of Ramesh et.al C113 to process fringes of any 
orientation. In this algorithm, fringe edges are detected first and 
any further detailed computation is done only within the fringe 
band. It is shown in reference 11 that a very simple method of 
global thresholding is sufficient to detect fringe edges. If there 
is intensity variation in the image, one can use either dynamic 
thresholdingC 1 1 3 or divide the image into convenient segments < tiles) 
and use appropriate global threshold for each of the segments! 11 3 to 
detect fringe edges. 

Once the edges are determined, the image is scanned row-wise <0 
deg. scan), diagonal-wise <45 deg. scan), column-wise (90 deg. scan) 
and cross— diagonal— wise (135 deg. scan). For each scan direction, 
the pixel having minimum intensity between the edges is selected as 
skeleton point. There may be situations wherein more than one pixel 
may have the same minimum intensity forming a band. If the band is 
of odd number of pixels, the center pixel is identified as the 
skeleton point and if it is of even number of pixels then two pixels 
at the center are identified. Thus at the end of initial processing, 
one gets four images of fringe skeletons corresponding to each scan 
direction. It^ is to be noted that these scans are done globally and 
are not for each individual pixel as it is done by Yatagai et.al. C83 
and Umezaki et.al.! 103. Noise points which are scan direction 
dependent, appear in each of the above scans. 

Further, the fringe skeleton is not continuous in one scan. The 
discontinuity is also scan direction dependent. 
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A novel approach with definite steps to retain which of the 
detected minimum intensity points as fringe points is reported by 
Ramesh and Pramod Cl 3 and the scheme is shown in Fig. 2. 2. 

!• The 0 deg. and 90 deg. scans are logically 'OR'ed to get a 
complete fringe skeleton. This image contains the noise due to 
both 0 deg. and 90 deg. scans. 

2. The 45 deg. and 135 deg. scans are logically 'OR"ed to get a 
complete fringe skeleton. This image contains the noise due to 
both 45 deg. and 135 deg. scans. 

3. Results of steps 1 and 2 are logically 'AND'ed to get the 
complete fringe skeleton free of noise. 

The procedure of identif ication of minimum intensity point in 
the above algorithm is developed with the premise that within a 
fringe, only one minimum intensity band can occur. The minimum 
intensity point is located directly without any intermediate storing 
of pixel values. However situations occur, particularly in the case 
of a saddle point, there can be more than one minimum intensity band 
along a scan direction within the fringe. In the improved version of 
the algorithm by Ramesh and Rajeev Singh C123 the gray level values 
along the scan direction within the fringe edges are first stored. 
The fringe skeleton points are then identified by using the minimum 
intensity criteria. By this approach, fringe skeleton points would 
not be missed from any number of bands. The above algorithm was 
originally implemented on the PIP system. One of the major tasks 
completed in the course of present work was to translate the program 
code from the PIP to the HVP-AT. The details are given in 
Appendix-l . 

2.4. Meed for a fringe tracing algorithm 

The fringe thinning algorithm based on the local minimum 
intensity point identification gives a fringe skeleton free of noise 
and discontinuities. The fringe skeleton has a thickness upto two 
pixels and has the intensity value of 250 i.e. (pure white). The 
skeleton is superimposed on the original image for clarity. 
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The fringe tracing algorithm has to perform the task of 
identifying pixels forming a particular fringe. The basic 
methodology is, if a fringe point is selected, the algorithm must be 
in a position to identify all contiguous points forming a fringe. To 
achieve this, an algorithm is proposed which basically uses a 3 X 3 
mask to search for the next fringe point. The intensity of the 
identified point is changed from 250 to some other value indicating 
the fringe order. For example 1 corresponds to fringe order 1, 2 
corresponds to the fringe order 2 and so on. The performance of the 
algorithm is tested for various class of test images. Following 
three classes of images were used in the present investigation. 

1. Simple geometric shapes generated using computer graphics. 

The methodology of image generation ensures continuous fringes 
of one pixel width. The images are selected to bring out the 
performance of the algorithm for all geometric properties of the 
fringes e.g. An ellipse is used to evaluate the performance for 
loops, the triangle and the rectangle are selected to test the 
performance for sharply bent fringes and a cross, to test the 
performance when the fringe has several branches a situation 
simulating a saddle point. 

2. Fringe skeletons extracted by the algorithm of Ramesh and 
Pramod Cl 3 on test images » 

The test image used for the purpose of present investigation is 
an image of concentric circles < ellipse - visually ). Figure. 2.3 
shows the fringe skeleton superimposed image. Though the final 
fringe skeleton is elliptic in nature, there is a difference between 
drawing an ellipse using the graphics utility and one extracted as a 
skeleton. The basic difference is brought out later in this chapter. 

3. Fringe skeleton extracted by the algorithm of Ramesh and 
Pramod Cl 3 on actual experimental fringe pattern. 




Fig. 2.3 Fringe skeleton extracted by algorithm o-f Ramesh 
and Pramod superimposed on the test image. 
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From analysing the results of the algorithm on the above class 
of images, the basic 3X3 mask used is improved upon. In the 
following paragraphs the details of the simple algorithm initially 
used and its improvements culminating in a complex fringe tracing 
algorithm is presented. 

2.5 Simple fringe tracing algorithm 

This algorithm uses a (3 X 3) mask as shown in Fig. 2. 4 The 
number represents the sequence in which the pixels are scanned from 
its local neighborhood for the fringe pixel. This mask has been 
selected giving priority to the pixels in the immediate 
neighborhood. The sequence of steps involved in tracing a fringe 
contour is summarised below. 

1. Starting from a point P(x,y> on a fringe <fringe order N) , the 
algorithm searches for the presence of a fringe skeleton pixel 
(having intensity value equal to 250) in its immediate neighborhood 
in the order shown in Fig. 2. 4. 

2. If such a point is located, then further search is stopped and 
the X and the y coordinates are set to the new values ( x',y'). 

3. The pixel intensity of the point P(x',y') is then changed from 
250 to the intensity value numerically equal to the order of the 
fringe which can be recognized by the program. 

4. The steps 1, 2 and 3 are repeated for the new pixel location 
P(x' ,y' ) . 

These steps are repeated till the end of the fringe is reached 
which is indicated by the absence of a fringe point among any of the 
neighborhood pixels within the mask. Having completed scanning along 
a particular direction, the X and Y coordinates are set to their 
initial value and scanning is repeated in the other direction thus 
tracing out the remaining part of the fringe. 
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1 6 



Fig. 2.4 Hierarchy of search directions for the 3X3 mask. 


2.5.1 Perfor»ar>ce for tl»e Images generated using computer graphics 

The mask performs satisfactorily for the rectangle. The mask is 
able to trace out the complete segments constituting the rectangle. 

In the case of the triangle, when the mask reaches near the 
vertex pixel it recognizes two pixels as lying on the fringe as 
shown in Fig. 2. 5a. 

The vertex pixel is skipped and the pixel on the next segment 
<direction 2) is selected as it has a priority over the direction 3, 
during the search operation Fig. 2, 5a. 

The mask is then advanced to this pixel. In the present 
position, fringe pixels exists in its immediate neighborhood one 
along direction 5 and the other along direction 8 <Fig.2.5b). The 
pixel along direction 5 alone is identified by the algorithm in view 
of the hierarchy of search used in the mask. Hence, the mask moves 
along that direction. In this position, all the neighboring fringe 
pixels are marked as shown in Fig. 2. 5c. The search routine 


terminates . 



Position of mask at 3 


Fig. 2.5 Performance of a 3 X 3 mask for image generated 
using computer graphics- 
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It can be seen that the mask is unable to trace out the whole 

fringe. If the mask is again setup at a point on the unmarked 
segment, then it traces out the remaining portion of the fringe. 
However , an approach like this makes the algorithm interactive 

rather than stand alone. 

In the case of a theoretically generated ellipse the algorithm 
gives a satisfactory performance. 

For a pair of intersecting lines, if initially a point is 
selected as shown in Fig. 2. 6a, the line is traced in accordance 
with the hierarchy set in the mask. On reaching the intersection 
point three fringe pixels exist along direction 6, 7 and 8 as shown 
in Fig. 2. 6b. The mask then moves only along direction 6 as shown in 

Lne Fig. 2. 6c. Thus, deviating from the original segment. The 

portion AOC is traced out and the portion BOD remains unmarked. 

2.5.2 Perf ormarice on fringe skeleton identified using the algorithm 
of Raraesh MTki Pra»c»d f 1 1 

The entire fringe is scanned and the overall performance of the 
algorithm is satisfactory. 

2.5.3 Performance for actual pl^toelastlc fringe pattern Images 

The algorithm does not perform satisfactorily and is unable to 
trace out some fringes completely. 

From the foregoing examples, it is seen that there are several 
situations wherein more than one fringe point pixel exists in the 
neighborhood. Ideally, all these fringe points need to be checked 
further. Some may be just noise points. If they are just noise 
points, a direction scan priority criteria may help in tracing the 
fringe. To achieve this, the direction of the fringe point for the 
previous pixel is to be stored and the existence of the next fringe 
point is to be checked along that direction first. 
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2. C Frlr»ge traclrig algorithm with direction scan 

1. The (3 X 3) mask is set up at the initial point and the next 
•fringe point location is determined and the point is marked. 

2. The direction of the fringe point from the current coordinates is 
then stored. 

3. The new coordinates are then set along this direction and the 
existence of a fringe point is then checked. 

4. If the fringe point is not found along that direction, only then 
the mask is setup and the next fringe point and the new direction 
are determined, otherwise the search is continued in the same 
direction . 

This algorithm works more efficiently in case of near linear 
fringes. The number of logical checks required, for such cases 
reduces from an average of four pixel point search to one pixel. 
Also, more priority is assigned for the previous scan direction 
rather than the priority for the positive Y direction which occurs 
in the mask. 

2.6.1 Perforeiance for the iHiages generated using computer graphics 

The algorithm performs well for the rectangle as well as the 
triangle. All the segments are traced out. 

In the case of the acute angle (Fig. 2. 7a) the horizontal 
segment is traced out till the vertex point is reached. Though the 
fringe pixel is also found in direction 2, since the previous scan 
direction is 3, it gets priority and mask traces out the complete 
segment in that direction. The fringe point is then located along 
the direction 8 and the tracing continues as indicated in Fig. 2. 7b. 
Thus the drawback of the simple fringe tracing algorithm is overcome 
with this modification. 



Position of mask at 3 


Fig. 2.7 Performance of a 3 X 3 mask with direction scan 
for image generated using computer graphics. 
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The algorithm performs well in the case of the theoretically 
generated ellipse. 

This algorithm gives a much better performance in cases of 
intersecting lines. As shown in Fig. 2. 8b, since priority is given 
to the previous direction ( direction 5 in this case), the correct 
part of the segment is traced out and as the mask does not deviate 
from the required path. However, the line segment COD will not be 
traced and it requires improvement to enable it to scan and trace 
all the branches. 

2.6.2 Perforinartee on fringe skeleton Identified using tl's? algorltl'im 
of Ramesh and Pramod 111 

The algorithm does not give a satisfactory performance with 
this image. The fringe is not traced out completely. This has been 
identified due to the problem of interlacing of pixels as shown in 
Fig. 2.9. 

A curved segment on the screen, observed under magnification 
shows the inclined segments like steps of a staircase. For a 
theoretically generated ellipse the best fit pixels are calculated 
corresponding to the position of the actual line segment. Hence they 
are represented by continuous (though staircase effect is evident), 
non interlaced pixels. At any point on the segment, only one fringe 
pixel exists in the neighborhood and the algorithm works 
satisfactorily. 

However, in the case of the actual fringe skeleton extracted, 
there occur more than one pixel representing the curve occur. This 
is because the pixel positions are located purely on the basis of 
the minimum intensity criteria. This leads to the mask straying in 
the wrong direction and causes the failure of the algorithm. As 
shown in Fig. 2. 9, the algorithm using the direction scan moves along 
direction 8 till the fringe pixel 4 . At this stage, the fringe 
pixels exist along directions 1 and 4, pixel 5. 




Position of mask at 3 


Fig. 2.8 Performance of a 3 X 3 mask with direction scan 
on a pair of intersecting lines. 
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Fig. 2. 9 


Per-formance of a 3 X 3 mask with direction scan 
for the case of the interlaced pixels. 

























In the view of the basic hierarchy of the search the mask moves 
along direction 1 and traces back along the second line of the 
pixels i.e. pixel 6 and thus deviates from the intended path. 

This problem of the mask straying from the expected path also 
occurs in the case of the actual photoelastic fringe pattern. 
However it performs better than the simple 3X3 mask. 

The next algorithm corrects the defect caused due to 

interlacing of the fringe pixels and proposes a more generalized 
approach for fringe tracing. 

2.7 Complex frli^ge tracing Algorithm t 

1. This algorithm also uses a 3X3 mask, with a difference that every 
pixel in the mask is scanned before a decision is made as to which 
direction the search should proceed. In this way all the branch 
points on the fringe are determined. 

2. If more than one fringe points are located within the mask, such 

points are stored to an array and the search continues in the 

direction as specified by the mask. Thus by the time the end point 
of the fringe is reached, all the branch nodes are stored in an 
array. 

3. When the end point of the fringe is encountered, the last point 
on the array is read, and the mask is set at this point and that 
part of the fringe is now traced. On encountering any additional 
branch points, they are stored to the array and the procedure is 
repeated till all those parts of the fringe are traced out. 

4. The count of the node points is decremented after the search is 

over. The mask is set at the previous branch point and the tracing 
continuous. Thus it can be seen that the branches are traced out in 

a reverse order i.e. the first node point being the last to be 

traced and so on. 



25 


a. 7.1 Periormartce for the images generated using con^^uter grapMcs 

In order to test the correctness of implementation, a tree 
image as shown in Fig. 2. 10. is generated using the graphics 
utility. The tree has explicit multiple branches from its nodes. The 
interlaced curve representation is thus a subset of this class of 
images and the performance of the algorithm is expected to be same 
for the tree image and the multiple branch points in an actual 
image. The node points are actually the ones which determine the 
next trace direction. 

The Fig. 2. 10 shows the performance of this algorithm for a tree 
structure. Starting from the lowermost branch, it reaches node 0 
first , stores the coordinates on the array, then continues along 
one path and reaches node point 1 whose coordinates are again 
stored. The scanning continues till the end of the fringe is 
reached. Now the mask is setup at node 1 now the node is removed 
from the memory and the other branch is traced till the node 2 is 
reached. This node is now stored as the last value and tracing is 
continued till the end of the fringe. The mask is now set to node 2. 
The node 2 is then removed from the array and remaining branch is 
traced out- This procedure is repeated for the remaining branches as 
they are called from the array and traced out one by one. 

The procedure also works well for the other types of 
theoretically generated images i.e. the ellipse , rectangle and the 
pair of intersecting lines. In the case of intersecting lines, the 
complete image is traced correctly. 


7- 2 Perf orniaince on fringe skeleton Identified using the algorltlini 
of Ramesh and Pramod £11 

The algorithm completely traces out the fringe skeleton in this 


case 




Fig. 2.10 Working o-f a complex -fringe tracing 
algorithm for a tree image. 


>(/) 
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The mections of the fringe having more than one pixel thickness 
(due to interlacing) are of same class as the fringe branches of the 
tree image. All such points are stored on the array. When the mask 
fails to find a fringe pixel ( it may be an end point or a 
discontinuity ) it returns to the last branch and continues tracing. 
The entire fringe is marked in this fashion. 

2.7.3. Perf ormaiice on fringe skeleton identified using the algorithm 
of Raroesh ar»d Pramod Cll 

The algorithm works well for the actual fringe pattern images. 
The mask returns to all the branch points and traces out the fringe 
points. A comparative performance of the previously mentioned 
algorithms for the image of a circular disk under diametral 
compression is shown in Fig-2.11. Though the complex fringe tracing 
algorithm is comparatively better, there are situations wherein the 
algorithm has not traced the complete fringe. On close scrutiny it 
was found to be due to fringe discontinuity at these locations as 
shown in Fig. 2. 12. 


2.8 Closure 

The relative performance of various algorithms show that the 
simple 3X3 mask alone is not sufficient for the actual fringe 
pattern images- The algorithm has also failed to trace completely 
the fringe even for some test images. 

The 3X3 mask with the directional scan, is computationally 
efficient and performs very well for the nearly straight line 
fringes. Complete fringes are traced out. However in the case of the 
curved fringes, if there is interlacing effect, the mask shows 
deviation and is unable to trace out the entire fringe. However the 
user can restart the mask by giving a point on the unmarked part of 
the fringe interactively which eventually makes the procedure 
semi-automatic . 
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Performance of the 3X3 mask for an actual Photoelastic 
pattern image. 



fringe 


Performance of the 3X3 mask with the direction scan for an 
Photoelastic fringe pattern image. 


actual 
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Performance of the complex fringe tracing algorithm for an 
Photoelastic fringe pattern image. 


actual 


Fig. 2.11 


Comparative performance of the 3 algorithms for 
actual photoelastic fringe pattern skeleton. 
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Fig . 



Magnified image of a discontinuity in the 
fringe pixels. 
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The complex fringe tracing algorithm works the best for actual 
fringe patterns. This algorithm is slightly less efficient than the 
directional scan algorithm but successfully traces out the complete 
fringes. This algorithm can be used for a wide variety of path 
tracing problems In actual implementation, the algorithm returns to 
all the branch points and scans for a fringe point even if that 
point is previously scanned by the other path. This drawback can be 
overcome by first determining if that branch is required to be 
scanned, after the end of the main branch is encountered. 

The algorithm stops at a discontinuity. In principle this 
drawback can be overcome by evolving a suitable criteria to check if 
no pixels are found at a particular stage whether it is actually a 
fringe end or a discontinuity. In view of lack of time it has not 
been possible to implement this feature. 

The next chapter discusses a new methodology for fringe 
ordering. The basic scheme, helps to order fringes along a 
particular scan direction. The procedure for fringe ordering 
discussed in the next chapter if coupled with fringe tracing 
discussed in this chapter can help to order the fringes in the 
entire field. 
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CHAPTER 3 

A NEW METHODOLOGY FOR AUTOMATING THE FRINGE ORDERING 


3.1 I ntrodijctl on 

Fringe ordering is one oi the crucial steps in the process of 
determining quantitative information from the fringe field. In 
general, the ordering of fringes has to be done systematically 
taking into consideration the physics of the problem and prior 
training is required to order the fringes. A few general guidelines 
to order the fringes are presented initially in this chapter. 

Automation of the ordering of Photoelastic fringes is one of 
the challenging tasks even today. Ordering of the fringes becomes 
easier if in the fringe field a zeroth order fringe exists and the 
fringe gradient directions are known. Ramesh and Pramod 113 
recognized the importance of the intensity data for fringe gradient 
calculation. They observed that at higher fringe orders, the grey 
level value of the minimum intensity pixel forming the fringe is 
higher. 

In this chapter, the observation made by Ramesh and Pramod 111 
on the intensity variation j.n the zone of stress concentration is 
further verified for a large number of images by plotting both 2d 
and 3d plots of the intensity variation. 

A new methodology is proposed to identify the fringe gradient 
direction using the intensity variation in the fringe field. Fully 
automated algorithm to order fringes in the case of a circular disk 
under diametral compression is developed. An interactive approach 
for processing crack tip fringe patterns is implemented and 
guidelines for processing a general class of fringe patterns is 
proposed towards the end « 



3.£ Fringe ordering in piiotoel Asticity ri3J 


Experience and intuition have made the numbering of fringes in 
common cases of experimental stress analysis, and in particular the 
numbering of isochromatics in photoelasticity, almost a trivial 
matter. In general, however, the problem is far from trivial. There 
is no systematic methodology for this purpose. 

Basically, most of the considerations can be linked to the fact 
that within the boundaries of a body the state of stress is uniquely 
defined and the stress and strains are continuous. Therefore on a 
line between two points with different values of stress; stresses 
with intermediate values must occur. Also, the stress tensor is 
uniquely defined so there can be only one state of stress at any 
point in the body. From this property, it follows that isochromatic 
lines of different order do not intersect. Because of continuity, 
the isochromatic lines are similar to the topological contours that 
end at the boundary of the pattern or form closed curves. The 
property of continuity helps in ordering the fringes. 

3.2.1 Guidelines to determine points at which shear Is zero 

When the white light is used, zero-order fringes appear black. 
This property does not help when monochromatic light source is used. 
Few of the geometric features which help in identifying the zeroth 
fringe order are presented below s 

Square Corners : 

Equilibrium requires that square outward corners on free 
boundaries have zero stresses in ail directions. They are therefore 
points of zero shear- Numbering may start from these points and 
increase continuously. Uniqueness of the stress field also ensures 
that in case of having more than one such fringe the fringe ordering 
can be carried out in any sequence and the fringes will have 
consistent order - 

Neighborhood of square corners s 

In many cases it has been found that the stress tangential to 
the free boundary reverses sign in the proximity of the corner. 
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The point at which this reversal occurs is therefore also a point of 
zero shear Fig.3.1. 

Reversal of sign on free boundaries : 

When there is enough knowledge of the boundary conditions and 
of the gross stress distribution in the body, to know that there is 
a reversal in sign along a free boundary, the point of change of 
sign must be a point where the shear is zero- The classic case of 
the circular hole in a plate subjected to uniaxial loading is an 
example . 

The Circular disk z 

The circular disk subjected to diametral compression is a 
unique case. All points on the circumference of the disk have zero 
shear . 

3-2^2 NLBsberlrig of fringes In the general case 

The uniqueness property of the stress values and the knowledge 
of one or several points at which shear is zero is sometimes not 
err‘ugh to determine unambiguously the order of all the fringes in 
the field. After giving certain fringe order it must be decided if 
the neighboring fringe order is one order higher or lower. 

Time derivatives : 

It is necessary sometimes to observe the development or the 

motion of the fringes, as the load is applied. As the load is 

increased, the first fringe appears at the point where r is a 

maximum in the field and then moves towards the region of lower 

T . When the applied load changes, all fringes except the zeroth 
max 

order one move. The point from which the fringes move away when the 
load increases, is at a higher fringe order than its surrounding and 
the point towards which the fringes move as the load increases, is 
at a lower fringe order. 

Space derivatives s 

An isDchromatic fringe pattern can be looked at as a 
topographic representation in three dimensions of the field of ; 
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The order of the fringes gives the height, or level of the maximum 

shear stress field with respect to the plane of reference. In 

Fig. 3. 2 certain points have been called "peaks** and "valleys** to 

help understand the analogy. The peaks are points where t is 

max 

maximum. Valleys are points where r is locally minimum. This 

max ^ 

minimum could be zero ( case of the 'isotropic' or singular points ) 
or could be a 'sink' where the stress is a minimum but has a finite 
value- In the case of the ring shown in Fig -3.2. there are four 
peaks on the inside boundary, at vertical and horizontal axis, and 
four valleys also on the inside boundaries, which happen to be of 
zero order. 

Inside the field of the ring, on the horizontal axis, there are 
two valleys (zero order) and on the vertical axis there are two 
other valleys at which the order is finite. 

Observation of the pattern in the ring shows that sometimes a 
fringe crosses itself or crosses one of the same order- This is 
obviously the case in the top half of the fringe pattern in Fig. 3-2. 
for the fringe of order 2.5. The order of this fringes can be called 
'transitional' order. The other transitional orders are those 

identified as about 1 and about 2 in the same figures. 

The intersection of two transition order takes place at a 
point, identified as S in the figure and sometimes called a saddle 
point- The transitional fringes crossing at point S divide the space 

in four regions- In two of these regions opposite to each other, the 

order of the fringes increases as the counting starts from S« In the 
other two regions the order of the fringes decreases. 

Saddle points are characteristically bounded by two 

isochromatic lines, one representing a shear somewhat greater than 
the one corresponding to the saddle point, and one with a shear 
somewhat smaller than the one corresponding to the saddle point. 
Both of these lines have two branches and the branches alternate 
with the branches of the other line to form a high-low-high-low 
pattern. Between the four branches there are asymptotic lines. The 
order of these lines was called 'transition order.' 



* 
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This property of saddle points permits one to decide whether the 
order to be assigned to isochromatic fringes should increase or 
decrease. 

3.3 Hew wethodology for automation of fringe ordering 

Automation of the ordering of Photoelastic fringes is one of 
the challenging tasks even today. Ordering of the fringes becomes 
easier if in the fringe field a zeroth order fringe exists and the 
fringe gradient directions are known. 

Identif ication of the fringe gradient direction is possible 
through any one or a combination of the following procedures. 

(i) Load the model gradually and study the fringes as they change 
with the load. 

<ii) Identification of a source and the sink in the fringe field. 
<iii> Use of white light which provides a colour code to identify 
the gradient as well as the fringe orders. 

From the point of view of automating fringe ordering, the 
procedure <iii> is a costlier option since the cost of a colour 
image processing system is at least three times the cost of a 
monochrome image processing system. The procedure <ii) is convenient 
for manual ordering of the fringes. However, if the information can 
be presented in a suitable form, it could be used for interactive 
ordering of fringes and in simple cases even complete automation is 
possible. Loading the model gradually is not the right choice for 
automation as one has to handle a large number of images to freeze 
the loading sequence. 

Changing the light source in two dimensional photoelasticity 
changes the fringe order at the point of interest and the model 
fringe constant is also correspondingly altered- Except the zeroth 
fringe order, other fringes shift when the wavelength is changed 
Umezaki et al C101 used this phenomena to identify the fringe 
gradient direction. The use of two different wavelengths has reduced 
the number of images to be processed to just two. 
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Ramesh and Pramod CID recognized the importance of the 
intensity data for fringe gradient calculation. They observed that 
at higher fringe orders, the grey level value of the minimum 
intensity pixel forming the fringe is higher. It is due to the fact 

that the monochromatic light sources do not have 1007. of the same 

o 

wavelengths. A monochromatic source of sodium vapour <5893 A) has 

c 

upto 4Z of another wavelength of 5896 A 11143. In the view of this, 
the ininimum intensity position forming the fringe is affected and 
the effect is pronounced at higher fringe orders. This leads to the 
umbrella shaped curve of the minima of the intensity values along a 
particular direction- Use of this phenomena requires only one image 
to be processed and is definitely advantageous. 

3.4 Use of two dimensioixal intensity plots for determining the 
stress gradient direction 

Figure 3.3 shows the plot of intensity variation along 
different lines in a circular disc under diametral compression. It 
can be seen that the variation in the minima of the intensity values 
is more pronounced at the higher stress gradient zones. This is 
because the monochromatic source used is not 100X of one wavelength 
as noted earlier. The value of the minimum intensity is higher as 
the point becomes more stressed. In a sense the model becomes more 
transparent at stress concentration zones. 

The same effect can be observed in the case of the crack image 
Fig. 3. 4 and Fig«3.5. 

The annular ring under diametral compression shown in Fig. 3.6 
provides a most general class of the fringes. This pattern has its 
unique characteristic fringes in the form of a source, the sink and 
the saddle points which have a positive fringe gradient in one 
direction and a negative fringe gradient in the perpendicular 
direction. If the ordering of the fringes can be automated for this 
case then one can say without exaggeration that any other fringe 
field can be automatically ordered. 
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Fig. 3.5 Variation in intensity for a crack tip fringe 
pattern . 
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Source and sink can be identified manually if the model is 
gradually loaded and the fringe movements are observed. Figure 3.6 
shows the intensity variation along different scan lines. The figure 
clearly shows that at higher orders the minimum intensity point 
forming the fringe skeleton has a higher value than that for the 
lower fringe orders. This observation again reinforces the fact that 
fringe gradient direction can be correctly identified in the 
neighborhood of the source. The question now is, can the intensity 
variation also indicate the location of a sink in the fringe field ? 
Figure 3.7 shows the variation in the intensity values along the 
line segment passing through a sink. It shows that the magnitude of 
the minimum intensity is lower for the inner fringe than that for 
its enclosing fringe, thus showing that the inner fringe is at a 
lower fringe order than the surrounding one, confirming the 
3 sink. This proves that the minimum intensity 
''^^i^tion can be used as a criteria for fringe ordering even in the 
neighborhood of a sink. 

3.5 Thi'ee dimei’islonal intensity plots for various types of fringe 
fields 

The two dimensional display is able to provide intensity 
information along a particular direction. However to obtain 
information of the whole field, a three dimensional plot of the 
intensity profile is drawn. 

Figure. 3. 8 shows the screen origin and the area of interest to 
be scanned for plotting intensity profile. 


Quadrant axis position 


Y 


^ X 




IMAGE 

WINDOW 


X 


y 


Fig. 3. 8 The image and the screen coordinate system. 
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Fig. 3.6 Variation in intensity along various scan lines 
for an annular ring under diametral compression. 



Fig. 3.7 


Magnified view of the intensity variation across 
a sink point in the annular ring. 
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Coordinate transforms are applied to a fringe pattern image 

with the X and Y values as the coordinates of the point < in pixels) 

and the corresponding pixel intensity as the Z coordinate. The 
image is scanned row wise. For each point the intensity value is 

read and assigned to the Z coordinate. The rotation transform is 

then performed and the image is displayed .The actual transform 
equations are as follows : 

1. The origin of the window is translated to the screen origin (T^, , 
T ). 

y 


X 


1 

Y 


i 


X - T 

X 

Y - T 

y 


2. Rotation about the X, Y and Z axis in that order, about the 

screen origin (R,R,R)is then performed. 

X y z 

R i.e. rotation about X axis by 6 . 

X X 


Y = Y cobB + 1 sine 

2 1 X 1 X 

Z = -Y sine + Z cose 

2 1 X 1 X 


R i.e. rotation about Y axis by 6 . 

y >' 


X cose - z sine 

2 y 2 y 

Y 

2 

X sine + Z_ cose 

2 y 2 y 



R i.e. rotation about Z axis by 6 . 
z 

X = X cose + Y sine 

4 5 S 5 

Y = -X sine + Y_ cose 

4 5 2 =2 
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3. The figure 


(-T 


K 



) . 


IB 


then translated back to the window origin 


X' = X + T 

4 X 

Y' - Y + T 

^ y 

4» The new X and Y values are plotted, giving a 3d projected view 
of the fringe pattern. 

In general, the rotation angles are so chosen that the final 
image obtained is aesthetically pleasing. The values ai B B and B 

X V T 

used in the present study was 10^^, -10"" and 45^^ respectively! Apart 
from these, a suitable scaling was also required in the Z axis 
coordinate. The values used in the present investigation range from 
0.65 to 1.75. Further , care needs to be taken while scanning the 
image. Finer scan was not beneficial and coarser scanning of the 
order of the 3 to 4 for normal as well as zoomed images, magnified 
by the order of 3 X 3 to 9 X 9 were found to be effective. 


Fig. 3. 9 shows the 3 dimensional plot of the intensity variation 
over the circular disk under diametral compression. Fig .3. 10 shows 
the zoomed view of the high stress gradient zone near the loading 
point. Figure 3.11 and Fig. 3. 12 show the intensity profile for 
crack-tip fringe patterns. 

Fig. 3. 13 shows the 3 dimensional intensity plot of the annular 
ring under diametral compression. Fig. 3. 14 and Fig.3.15 shows the 
zoomed view of the sink point and the saddle point. The 
topographical features help in easy identification of the variation 
in the fringe order. 

The 3D plots have clearly brought out the fact, that there is 
substantial intensity variations in the fringe field. Use of 
intensity information has not been paid much attention earlier, in 
view of the difficulties in recording intensity information. 
Intensity was monitored at a few selected points or along specified 
lines using a photo multiplier tube. 















3.11 3— dimensional inten 
tip fringe pattern 




12 3-dimensional intensity profile for 
tip fringe pattern with a frontal loop 





Fig- 3.13 3-dimensional intensity profile for annular ring 
under diametral compression. 


Fig. 3.14 3-dimen3ional plot of the magnified view of the 
sink point. 
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Fig. 



saddle point. 
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Collecting intensity data over a field was unthinkable earlier. The 
modern developments in the technology have given rise to the 
development of CCD cameras, which can quantize the intensity values 
over the entire field at video rate. The resolution of the cameras 
are substantially higher to capture even minute variations in the 
intensity information. This development has showed promise to use 
intensity variations for data processing. 

3.6 Automation of fringe ordering for special class of fringe 
patterns 

It has been established in sections 3.4 and 3.5 that intensity 
variation in the fringe field could be used to judiciously identify 
fringe gradient direction. To achieve this, a proper choice of scan 
direction need to be selected. The problem becomes much simpler, if 
the nature of the fringe field is known apriori. 

Knowing the gradient direction, if fringe order of any point 
lying in the gradient direction is known, then it is possible to 
identify fringe orders for any other point. To know the ordering of 
the fringe field, fringe tracing algorithm must be used to trace out 
the complete fringes. The new scheme involves the following steps. 

1. Fringe thinning to identify fringe skeleton. 

2. Fringe ordering along a particular scan direction using, 

Ci) User specified basic informtion on the model geometry. 

(ii) Intensity variation over the field to identify the fringe 
gradient direction either independently or along with step (i). 

3. Use of fringe tracing algorithm to order the complete fringe 
field. 

The scan path, has to be found out such that all the fringes 
within the field are located and scanned in proper order. This fact 
is very critical to ensure correct ordering of the fringes. This 
information has to be provided by the user in case of the general 
class of the images. Knowledge of the model geometry can provide the 
following information. Location of the crack tip in case of the 
crack problem, or coordinates of the loading point enables the 
program to determine the "source". 
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If the free Burf.ce Ib indicated then gradient direction can be 
established. For specific class of problems, i.e. the circular disk 
under diametral compression and the mixed mode fringe pattern, the 
scan direction have been determined to work for most of the fringe 
patterns of these classes. 

3.7. Automation of fringe ordering for a circular disc under 
diametral compression 

Calibration of photoelastic materials is one of the important 
step in any photoelastic experiment. Many laboratories use circular 
disk under diametral compression as a model to determine this. 
Automation of fringe ordering and subsequent analysis of data for 
such a case is desirable. 

Circular disk comes under a special category as noted in Sec. 
3.2.1. that the boundary of the disk except the loading points is at 
zeroth fringe order. The fringe pattern is symmetric about both X 
and Y axes. The fringe order increases from the free surface, to the 
loading point. Hence the scan path must, ideally be perpendicular to 
every fringe and should start from the free surface and go upto the 
load point Fig. 3. 16. The positive X and V axis are a correct choice 
for the scan paths. The algorithm searches for the fringes from the 
farthest fringe point on the positive X axis to the centre of the 
disk and from the centre of the disk, along the positive Y axis till 
the load point. It must be ensured that the loading points should be 
more or less oriented along the Y axis. Implementation of the 
algorithm requires the determination of the centre and the diameter 
of the disk. 

An interactive procedure is implemented to identify the centre 
and the boundary of the circular disk and a circle is drawn using 
the graphics command, corresponding to this centre and radius 
chosen. If the circle lies exactly on the boundary of the disk then 
these centre and radius values are accepted, otherwise the circle is 
erased and the procedure is repeated. 

1. The image is cleaned and then skeletonised using the algorithm by 
Ramesh and Pramod C13. This image is loaded in quadrant 0. 
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2. Since the boundary of the disk is known, the farthest fringe 

point on the X axiB is known. 

3. The i and the j values are initialized to these coordinates and 
the complex fringe tracing algorithm is used to trace this fringe. 
The fringe is marked with intensity value equal to its fringe order. 

4. On reaching the end point of the fringe, the original i and the j 
values are read and the pixels are read along the scan path to 
locate the other fringes. 

5. The mask is moved along the X axis, towards the centre of the 
disk till the next fringe is located. The fringe order is then 
incremented and this fringe is traced out using the complex fringe 
tracing algorithm . 

Steps 4 and 5 are repeated till the mask reaches the centre of the 
disk . 

6. The previous fringe order is then stored in the memory and the 
mask continues to scan for the next fringe along the negative X 
direction, away from the centre. The fringe orders are now known as 
the fringes are also symmetric about the Y axis. 

7. When a fringe point is located, the fringe is traced out in the 
intensity corresponding to its fringe order and the fringe order is 
decremented . 

This procedure is repeated till the disk boundary is 
encountered. Now the mask is set at the disk centre and the other 
part of the scan path is searched. 

8. The fringe order stored in the memory is recalled and incremented 
by one. The mask now scans for a fringe point along the positive Y 
direction. The fringes are marked as the fringe points are located 
and the fringe value is incremented by one. This procedure 
is repeated till the the boundary of the disk is reached. 
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9. Same procedure is repeated -for scanning the fringes along the 
negative Y direction upto the other loading point, thus ordering all 
the fringes. 


10. The fringe ordering obtained by the above procedure is cross 
checked by comparing the intensity profile gradient to fringe 
gradient - 

3.8. Interactive fringe ordering for the mixed mode crack-tip fringe 
pattern 

Selection of the scan path : 

The mixed mode crack fringe patterns are non symmetric. Hence a 
more complex criteria is required to determine the scan path. The 
problem is split into two parts determined by the geometry of the 
fringes. The fringes are in the form of either elliptical loops 
oriented along the Y direction or near straight fringes ( parts of 
big loops which go out of the image field ). Apart from these, in 
some cases frontal loops also exist. 

The scan paths are different for each of these classes of 
fringes. For the elliptical fringes oriented along the Y direction 
and the linear near vertical fringes , a line parallel to the X axis 
(the crack axis) and at sufficient distance from the origin on 
either side is selected, as indicated in Fig. a. 17. The distance is 
decided by the resolution of the fringes close to the crack tip. 
This value is usually about 2 to 3 times the crack thickness as seen 
on the screen i.e. 10 to 15 pixels. The line segment is from the 
negative X extent of the image window upto the crack tip. The scan 
path is required to extend beyond the crack tip to trace the open 
loops or the near linear vertical fringes on the positive X axis 
(Fig. 3. 17). For the closed and the open frontal loop , the scan path 
is selected along the Y direction at some distance on the positive X 
axis, from the crack tip. This scan path is from the positive extent 
of the Y axis to the crack tip Y coordinate. 
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In general » one does not have a zeroth fringB order along these 
directions, hence the user is required to provide the -fringe order 
for one of the points on these fringes. 

Algorithm t 

1. The image is cleaned and then skeletonised using the algorithm by 
Ramesh and Pramod Cl]. The original image is loaded in quadrant 0. 
It is to be ensured that the orientation of the crack should be 
along the positive X axis. The crack tip is to be located by the 
user, with the help of the cursor. 

Scan Path 1 : 

2. The scan path is set up from coordinate corresponding to the 
negative X extent and Y coordinate is 10 pixels above the crack tip. 

3. The cursor is displayed and the user is asked to input the 
fringe order of the first fringe intersecting the scan path. 

4. When a fringe point is located, the complex fringe tracing 
algorithm (section 2.4.4) is used to trace out that fringe then the 
fringe order is incremented and the scanning continues till the 
crack tip X coordinate is reached. 

Scan path 2 : 

5. The knowledge that ahead of the crack tip along the crack axis a 
sink is present, is used. The second scan path, parallel to the Y 
axis is set and the user is asked to input the fringe order of the 
outermost fringe along this path. The pixels along this path are then 
scanned for the frontal loop and the open fringes. 

6. On reaching the fringe pixel the fringe is traced out and the 
value is decremented by one. This is because of the presence of the 
sink . 


This procedure is repeated till the crack tip Y coordinate is 
reached. 
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Scan path 3 i 


7. Now the open loops and the near linear fringes are traced. The 
pixels are scanned along the path parallel to the X axis, from the 
crack tip to the X window extent. 

8. When a previously traced pixel is located the fringe order is set 
to this value. 

9- On encountering a fringe pixel the fringe order is decremented by 
one and it is traced out. 

Thus the complete set of fringes around a crack tip are scanned. 

In view of the interactive selection of scan paths the procedure is 
only a semi-automatic. 

3.9. Closure 

It has been successfully demonstrated in this chapter that 
intensity variation in the fringe field is not trivial but can yield 
a very useful information on the fringe gradient directions. Though, 
fringe ordering for complex image is difficult, a few general 
guidelines to order them are presented. 

While developing the automatic fringe ordering algorithms for 
circular disk and crack tip images only minimal information need to 
be input by the user. Further, the information needed was simple, 
and any one with basic background in Mechanical Engineering will be 
able to provide this information. 

The annular ring image represents one of the most complex 
fringe pattern. This image is used as a further check to utilize the 
concept of using intensity variation for fringe gradient 
identification. Typical intensity variation in the neighborhood of a 
'‘source** and a "sink” were obtained. Following the lines of the 
circular disc and the crack-tip image it is possible to develop a 
module to identify fringe ordering in an annular ring. 



57 


FrinQP ordering is only an intermediate step and one needs 
methodologies to extract field parameters using the fringe data. The 
next chapter provides a methodology to extract the stress field 
parameters in the neighborhood of a crack tip from the information 
of the fringe order and positional coordinates of various points. 
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CHAPTER 4 

EVALUATION OF STRESS INTENSITY FACTOR USING A HIGHER ORDER TICC^Y 


4. 1 Intr oducil on 

One of the important stress field parameters that is to be 

evaluated in Fracture Mechanics problems is the Stress Intensity 

Factor(SIF). Photoelasticians recognised way back in 1958 itselfClSl 

the fact that atleast two parameters are needed to characterise the 

stress field in the close vicinity of the crack-tip. The correct 

determination of is possible only if both the parameters namely 

and are determined from the stress field. The ingenuity of 

Irwin C153 to generate two equations from a single data point from 

the fringe field, made the evaluation of Mode I SIF and a very 

cx 

simple and did not warrant the use of digital computers. 

Sanford C33 proposed an overdeterministic method to evaluate 
and cf in a least square sense using information from the field 

wX 

rather than a single point as was used by IrwinClSl. The data to be 
collected included the fringe order and positional co-ordinates of 
various points in the field. In view of the large volume of data to 
be collected, there was increasing demand at this stage to automate 
the data collection process. 

Most of the practical structures of engineering importance have 
finite dimensions and the cracks in them may be located in zones of 
stress concentration and the cracks may be large enough so that the 
crack-tip may be closer to a boundary. The basic stress field 
equations of fracture mechanics are developed for cracks lying in a 
body of infinite dimensions <i.e., the crack-tips are far away from 
the boundaries) subjected to an uniform stress field. When these 
equations are attempted to be applied for solving . eai life 
problems, one gets erroneous results- 
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The recognition of this tact has lead the researchers to 
account tor this ettect and mul ti— parameter stress tield equations 
have been developed in the past two decades. 

Among the muti-parameter stress tield equations, the simplicity 
equations proposed by Atluri and KobayashiC21 are the simplest one 
and their equivalence with other equations is brought out in Ret. 11 
Using these equations and tol lowing the methodology ot SantordESl, 
an overdeterministic least square method is developed to extract the 
stress tield parameters trom the isochromatic tringes. 

4.2. Developments In the description of the stress field equations 
in the nel gidior hood of a crack-tip 

4.2.1 Single Parameter Stress Field Equation 

The Westergaard C163 complex stress tunction technique tor 
solving opening mode crack problems has played an important role in 
the development ot linear elastic tracture mechanics. Westergaard in 
his paper gave the stress tunctiohs tor the cases ot a central crack 
and a series ot equally spaced straight cracks ot length 2a in an 
intinite plate with a biaxial tield ot tension a. Westergaard 
originally reported that his stress tunction models the stress tield 
tor the case ot a crack/trracks in an uniaxial stress tield. However, 
it has been shown by later investigators that Westergaard ' s solution 
is valid only tor a biaxial tield. IrwinClTl in 1957 tol lowing the 
semi— inverse procedure suggested by Westergaard added three 
additional examples. They are 

(1) Single crack along the x-axis extending trom -a to a 
with a wedge action applied to produce a pair ot "splitting torce 

ot magnitude P located at x = b. 

(2) The situation ot above example with an additional pair 

ot torces ot magnitudes P at x = -b. 

<3) Exampled) repeated along the x-axis at intervals Z 

and with the wedge action centered so that b is zero. 
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IrwinC173 showed that in all the above five cases the stress 
distribution near the end of the crack can be expressed independent 
of the type of loading if the region of interest is very close to 
the crack— tip (i.e. r/a and r/(a— b) may be neglected in comparison 
to unity) and introduced the concept of Stress Intensity 
Factor(SIF). The expressions thus obtained are known as classical 
Westergaard equations and for a Mode I loading situation they are 
given as, 
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d. 2. 2. Two Parameter Stress Field Equation 


Irwin C173 observed that when an extending crack moves across a 
plate of finite width, the crack may attain sufficient length so 
that the stress— field equations obtained from the consideration of a 
crack in an infinite body may not be accurate. He reported that 
correction needs to be applied for a finite body and introduced a 

term. Thus, the modified 
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In a discussion of the paper by Wells and PostC183, IrwinC153 
in 1958 reported a method to evaluate, as well as stress term 
from the isochromatic fringe patterns. Using these values he 
theoretically plotted the contours of maximum shear stress 
(isochromatic fringes) and showed that the prominent tilt seen in 
the experiment is well represented. 
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On the otherhand if one uses classical Westergaard ' s equations one 
gets the fringes symmetric with respect to both x and y axis. This 
situation occurs in the experiment in a region very close to the 
crack~tip. In a sense the use of the second term has helped in 
enlarging the zone of data collection. 

Sih C193 in 1966 starting from the Goursat— Kolosov complex 
representation of the plane problem showed that the symmetry 

condition t = 0 on y = 0 could be satisfied by a less restrictive 

xy 

assumption than employed by Westergaard. He obtained near — field 
stress equations which included a constant term to a term. Later 

X 

Eftis and Liebowitz C203 in 1972 showed that Westergaard ' s stress 
function can be related to the Goursat— Kolosov complex 
representation. They confirmed that Sih's equations are the same as 
Irwin's modification to the Westergaard equations. 

4.2.3. Multi Parameter Stress Field Equations 

Sanford in 1979 C33 observed for the first time that when the 
crack length increases in an Modified Compact Tension (MCT) 
specimen, several isochromatics cross the crack-axis <Fig.4.1). 
Sanford reported that this phenomenon is due to the finite 
dimensions of the specimen. In order to account for finite bodies, 
Tada, Paris and Irwin C213 introduced a series form of Westergaard 
stress function. Sanford C223 observed that the modified Westergaard 
functions predict only a constant fringe order along the crack— axis. 
However, for cracks approaching a boundary, several fringe orders 
cross the crack- axis. This counter-example raises questions about 
the validity of the modified Westergaard equations. 

Starting from this standpoint, Sanford reported that an 
additional stress function Y(z) should be added to the series form 
of Westergaard stress function Z(z) where z is x+iy and i is -/=!. He 
showed that the condition t = 0 on y = 0 can be achieved by a less 
restrictive assumption than used by Westergaard. 




Tinge pattern observed 
specimen. Isochromatics 
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He selected the additional stress function Y(z) such that it 
satisfies the condition Im(Y( 2 )) «= 0 on y = 0, He expressed the 
stress function Y(z) in a series form. With this additional stress 
function, Sanford and his co-workersC223 showed that crack-tip 
stress fields can be better characterized. In all the papers 
reported by Sanford and his CD-WDrkersC223 the stress field is 
represented only in terms of complex stress functions and they are, 
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where Z' and Y' are the first derivatives of Z and Y 
respectively with respect to z. 


For a single-ended crack, with the origin of co-ordinates at 
the crack tip and the negative x— axis coinciding with the crack 
faces, the functions Z(z> and Y(z) can be expressed as; 
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where C and C are real coefficients. 

2; 2jfl 


(4.4) 
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(Note : In the literature series coefficients of Z(z) are denoted fay 

A. and for Y(z) as However replacing them by ^ and 

ii simpler to collect powers of r in the series solution of 


stresses) 
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One can obtain stress field equations in terms of r ,& by 
substituting Eq.<4.4) and Eq.(4.5) in Eq.(4.3) and expressing z in 
terms of r and t). However, the resulting equations are lengthy and 
can be found in reference 22. 

Equation (4.3) is known as the Generalised Westergaard 
equation. The classical Westergaard equations are obtained by making 
Y(z) zero and Irwin's modified equations are obtained if the 
function Y(z) is taken to be a real constant. 

Atluri and Kobayashi C2] in 1987 have reported the stress-field 
in series form explicitly. The stress field equations reported by 
them are compact and is elegant to use for devising least square 
algorithms for evaluating stress field parameters. The stress field 
for the mode I case is described as follows: 


cos(J-3)e‘ 
cos(^-3)© ^ 
sin (^3)0 

(4.6) 

where A = K / Y2nr . (Note: In the original reference it is reported 
11 I' 

as K^y2frr which is a typographical error) and 4 A^^ = 

These equations have been applied to a practical mode I problem 
to test their validity and is discussed latter. 
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4.5 Mixed Mode stress field equations : 

Atluri and KobayashiC23 also reported the stress-field 
equations for the mixed mode case. They are described as follows : 
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where = Kj VSnF and (Note; In the original 

reference there was a typographical error in the term of the mode 
II series. The first term was a 'sin' term and the second term was 
denoted as a cos term. It was found in the course of present 
investigation that the second term must also be a 'sin' term. This 
was made possible with the help of graphics software which showed 
symmetry in fringes only when the term was changed to sin . It was 
also counter checked by comparing the first term of the series 
solution to the singular solution reported in the literature.) 

4.6. Developments In SIF evaluation methodology 

The Irwin's method for SIF extraction from photoelastic 
patterns was the accepted method for analysis for many years. In 
this method the positional co-ordinates and fringe order 
specific point that satisfies the criteria /d© = 0 are used for 

evaluating Kj and Kobayashi and Bradley C231 modified the 

formulation, but fundamentally, the approach was unaltered, 
then several investigators reported various methods and Etheridge 
and Dally C24] reported a review on the relative accuracy o 
various two-parameter methods. The review brought out the fact that 
the zone of data collection has to be very small of the order of 


r/a < 0.03. 
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Though photoelaetici ty gives whole field information, the 
previous reluctance to use information from the entire field stemmed 
primarily from practical limitations on data reduction. The 
positional co-ordinates of each data point used had to be measured 
manually and the results are calculated with the aid of a hand-held 
calculator or with a main-frame computer. This approach generally 
required manual data entry. With the introduction of small 
microcomputers and digitising tablets and the more recent PC based 
image processing, the drudgery of data reduction is removed C253. 

Improvements in data collection methodology have lead to the 
development of several techniques to process the data from the field 
rather than restricting the analysis to a point (as in Irwin's 
method) or a few points along a specified line. Among the several 
methods C263 the ovei — deterministic approach of Sanford and Dally 
[273 has come to stay as the most useful technique for evaluating 
the SIF and o- . 

OX 

All the approaches described above share a common feature, 
namely that the analysis procedure used to extract the SIF from the 
full field fringe patterns relied on the near field equations of 
fracture mechanics. In principle, the accuracy of these methods 
should improve as the region of data acquisition is reduced to 
smaller regions around the crack-tip. However, it is generally 
recognised and has been independently verified by several 
researchers 128,293 that the plane-stress assumption ceases to be 
valid in a small region at the crack-tip. In this region the stress 
field is three dimensional and experimental observations are 
influenced by stress gradients through the thickness. In addition, 
the stress field is altered by localized crack-tip blunting. On the 
other hand, if this region is excluded from the analysis, the region 
of data acquisition may not lie fully within the domain in which the 
inverse— square root singularity dominates the fringe pattern 
(Fig. 4. 2) . 




Fig. 4.2 Singular non singular and far field region around a crack 
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The size and shape of the singularity dominated zone is an important 
factor to be determined and such studieBC30,31 3 emphasized the fact 
that the singularity dominated zone is very small and indicated the 
use of a multi-parameter stress field equations for data reduction. 
The zone of data collection in such a case can be larger and one can 
collect data from the mixed field zone also (Fig. 4. 2). 

4.7. Evaluatlori oi stress Iritensity factors using least square 
teclinlqoe 

Stress-optic law relates the fringe order N and the principal stress 
as 

N F 

= i a - <y> ( 4 . 8 ) 

t ^ ^ 

where is the material fringe value and t is the model thickness. 
For a plane stress problem the principal stresses are. 
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Substituting Eq.<9) in Eq.(8) , a function g is defined for the 

data point as follows 
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If Eq.(4.7) is substituted in Eq.(4.10), then the Eq.<4.10> is 


non-linear in terms of the unknown 
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where the subscript 'i' refers to the i^*^ iteration step and 

• '“rn'i’ ‘“ija’i'' •'"'iii’i 

denote the corrections to the previous estimates of ^^I2^i’ 

• ■ ’ ’^^iii^’ ^^ii2^i’ ■ ■ 


Corrections are determined such that ^^m^i-M” ® thus 


Eq.di) gives 
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Applying this iteration scheme to M ( M > k + / ) data points 
results in an overdetermined set of linear equations in terms of the 

unknown corrections ^^I2^i’ * ’ ’ ^^ik^i ^^iii^i’ 

(AA ) , ... , (AA ,) . given in the matrix form by, 

II2 t’ lit t 
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(4.14) 
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In forming the Eg. (4.13) one needs to calculate 

(dg^/dAj^) , (dg^/dA^^) ... etc. and at first it appears that it may 

be difficult to do it. However , one can determine them very easily 

as follows. A typical derivative for the data point with respect 

to A is obtained by taking the derivative of Eq.(4.11). 
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The Convergence can be decided by either of the two following 
criteria, namely (1) the parameter error minimization and (2) the 

fringe order error fwinimization . 


In the parameter error minimization the iterations are stopped when 
the values of •• » 

reasonably small ( say of the order 

of 10~^) - In the fringe order error minimization, using the newly 
calculated values of 

-fringe orders corresponding to the 
selected data points are calculated theoretically during every 
iteration step and are compared with the experimentally obtained 
fringe orders. The convergence criteria is satisfied if , 
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The convergence error is usually of the order of 0.05 to 0.1 


A software has been implemented on the PC-“3S6 environment in 
Turbo Pascal for the least square analysis and the theoretical 
reconstruction using upto 13 terms each of model and mode II of the 
Atluri and KobayashiC23 form of the stress field equations. The 
software requires the data file containing the X,Y coordinates and 
the fringe order. Converts the data in real world coordinates, using 
the data preprocessor, then the least square procedure is invoked 
and the results are written to another data file. Using these 
parameters the program recalculates ^ 2 ^ values and hence the 
fringe order for every point in the data field. The fringes are then 
plotted using the graphics utilities of the Pascal environment. 

4.8. Performance for a mode I crack problem 

Figure 4.3 shows the isochromatics observed in an annular ring 
subjected to internal pressure with cracks emanating radially from 
the outer boundary. 
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Fig. 4.3 Isochromatic fringe pattern obtained for an annular ring 
subjected to internal pressure with cracks emnating from 


the outer boundary. 
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The fringe patternB observed close to the crack-tip are similar to 
the one observed in MCT specimen for large crack lengths. This 
example shows that this feature is due to a combination of two 
aspects namely the finite dimensions of the specimen and also due to 
the fact that the crack is lying in a zone of stress gradient. 

A least square analysis is carried out to determine the minimum 
number of terms required to model the stress field accurately. 

The minimum number of data points required for the least square 
algorithm to work satisfactorily [321 is given by 

M = 2(k + /) + 2 

where k and C are the number of Mode I and Mode II terms and M is 
the number of data points. Sixty five data points were collected 
from the sketetonized fringe pattern using the inhouse developed 
software. The points were selected so as to represent the full 
characteristics of the data field, i.e. Depending on the size of the 
fringe, equally spaced data points were collected from the field. 
The points selected are also skewed with reference to the X axis. 
This ensures that the least square equations do not become ill 
conditioned . 

The convergence criteria selected is the fringe order error 
minimization. The parameter values are resubstituted in the 
equations and the isochromatic fringe pattern is theoretically 
reconstructed. This is compared with the actual fringe pattern by 
echoing back the data points collected for the analysis. The 
analysis is performed using upto 8 terms of the series solution and 
the results are as shown in Fig. 4. 4. 

Figure 4.5. shows the graph of the cumulative percentage error in 
fringe order with the number of terms of the least square solution. 

The percentage error for the 2 term and the ^ term solution is 
about 24%. 
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Fig . 4 



5 Graph of cummulative percentage error v/s number of terms 
of the series solution. 
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The parameter values are not able to simulate the general stress 
field. This is obvious, since the first 2 terms are valid only for 
the singularity region and not for the general stress field. 

With the inclusion of the fourth and the fifth term, the 
accuracy of the solution improves drastically. The total percentage 
error in the fringe order comes down to about 67.. The series 
solution is able to model the general stress field with a better 
accuracy, however it is still not able to model the inner frontal 
loop. 


The solution stabilizes for 6 or more terms in the series 
solution. The percentage error in the fringe order drops to about 2 
percentage. The stress field is modelled accurately. Beyond 6 terms 
the accuracy improves only slightly. The 8 term solution is the best 
one for this class of fringe pattern. The computational effort 
required for a least square solution increases in the order of the 
square of the number of parameters used. Hence the computational 
time required increases drastically as the number of terms 
increases . 

Figure 4.6. shows that the KI value stabilizes with the 
increasing number of terms used. Figure 4.6 also shows the variation 
of the other coefficients of the stress field equations as a 
function of the number of parameters. 

4.9. Closure' : 

In this chapter the use of the equations proposed by Atluri and 
Kobayashi C23 for evaluating the stress field parameters in Fracture 
Mechanics problems is discussed. An overdeterministic method for 
evaluating stress field parameters in a mixed mode situation is then 
presented. The result obtained shows that the program is validated 
^or solving a Mode I situation. In the next chapter the use of the 
above methodology for evaluating stress field parameters for a 
mixed-mode situation is discussed. 
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CHAPTER 5 

EVALUATION OF STRESS FIELD PARAMETERS IN THE NEIGHBORHOOD OF A CRACK 
LOCATED IN THE TENSILE FILLET OF A SPUR GEAR. 

5.1 Intr oductlon 

Gear is one of the most important and most widely used machine 
element. The Qear tooth failure commonly occurs due to fracture 
initiated due to a fatigue crack at the tooth fillet. There could be 
several microscopic flaws during manufacturing or developed during the 
actual service conditions. They grow in size during the operation and 
can reach critical size, leading to failure. The effect of the flaw on 
the overall strength of the machine element has to be evaluated. If the 
stress field around the crack is known then using the principles of 
fracture mechanics one can evaluate the residual life of the gear or 
the reduced load to prevent crack propagation. 

There has been only a limited study in the literature in analysing 
the stress field in the neighborhood of a crack at a tensile fillet of 
a spur gear. Honda and Conway C333 attempted to study this problem 
using the Finite Element Technique. They modelled only one gear tooth, 
assuming that it is rigidly fixed to the gear hub. They observed that 
the stress field is a combination of both mode I and mode II loadings- 
They evaluated the S.I.F for various positions of load application 
points from the addendum side to the dedendum side. 

There has been no experimental verification of the results 
reported by Honda and Conway C333. They basically analysed only one 
gear tooth since modelling the entire gear is difficult from the memory 
requirement point of view unless the concept such as cyclic symmetry is 
invoked. The analysis can be considered realistic only if it 
incorporates the loading due to assembly and the interaction of 
adjacent gear teeth. Photoelastic analysis has certain advantages in 
terms of the accuracy in simulation of the actual loads occurring on 
the model, like the clamping load, the hub stiffness etc. 
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In the previous chapter the use of the equations proposed by 
Atluri and Kobayashi C23 for evaluating the stress field para.-.-L— s 
in Fracture Mechanics problems were discussed. In this chapter the 
application of these equations to the mixed mode crack problem is 

presented . 

5.2: Experl details 

Photograph in Fig^S-lpa shows a view of the loading fixture as 
mounted on the loading frame of the polariscope during the experiment. 
Figure S.l.b shows the closeup view of the loading arrangement- The 
gear loading frame comprises of two shafts with a variable center to 
center distance. The top shaft is set by means of the protractor scale 
attached on it and rigidly clamped at the required position. The 
alignment of the shafts is adjusted by means of the spherical washers 
and the locating strip C343. The prime objective of the alignment is to 
maintain the loaded gears in the same plane every time they are loaded. 
When the axes of the two gears are exactly ' paral lei , then the locating 
strip moves freely in the guide block. The gears are keyed to the shaft 
and clamped by means of lock nuts. The bottom shaft is free to rotate 
in its housing. To the other end of the bottom shaft is keyed a 
pul ley. A wire rope is wound over it and to the other end of the wire 
rope is attached a pan. Suitable load is applied using this 
arrangement. 

A pair of spur gears made from a precast sheet of Araldite sheet of 
thickness 6.26 mm having the following properties is used : 

Modulus of Elasticity, E = 2.6 SPa 

Poisson ratio, !> = 0.33 

Spur gears are cut with the following specifications » 

Pressure Angle = 14.5 

Diametral Pitch = i> 
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Fig 


1 .b Closeup view 


o-f the gear loading -frame 
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(Module 
Gear Ratio 
Number of teeth 
Axial face width 


- 4.3333 mm) 
>= 1 
= 24 

= 6 . 26 mm 


5.3 Selerctlon of lf»e crack location 

In the S.I.F. analysis of the crack at the tensile fillet of the 
gear tooth using Finite Element Method, Honda and ConwayC333 reported 
that the crack tip had been selected at the highest stressed point in 
the tensile fillet region of the gear tooth. To identify this point 
they obtained the stress field in the fillet area initially. Following 
a similar approach isochromatics are recorded for the gear tooth fillet 
using photoelasticity. It can be seen clearly form photograph in Fig. 
5.2 that the isochromatics faithfully follow the fillet contour. In 
other words, the value of the maximum shear stress is almost a constant 
along the entire fillet. Thus there is no specific point along the 
fillet, where there is maximum stress. In a sense spur gear is designed 
with a contoured fillet which is good from service point of view. The 
mid point of the fillet is chosen as a probable location for a crack. A 
thin slot of 1.26 mm long, and 0.1 mm thick is cut simulating a crack 
at this point normal to the surface. 

5.4 Frli’ige ordering 

The fringe ordering is a difficult procedure in this case because 
of the complexity of the stress field. The time derivative of the load 
and the colour code are used as criteria for the fringe ordering 
(section 3.2.2). The fringe patterns are recorded and the direction of 
the motion of the fringes is studied as the load is increased to 
determine the fringe ordering. Figure 5.3 shows the photograph of the 
fringe patterns for different loads under monochromatic and white light 


sources. 
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Fig. 5.2 iBDchromatic Fringe pattern o-f a gear tooth under 
under monochromatic and white light. 
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The colour code Bpecifies that the fringe that is black is a 
zeroth order fringe. Fringe orders 1, 2 etc are identified as red-blue 
and red-green transitions. The second order is identified by pinkish 
red and the third order is identified by orange red. These colours are 
easily identifiable when the isochromatic fringe field is recorded 
using white light. The corresponding pictures taken using a 
monochromatic source ( Sodium Vapour ) shows the fringes as black along 
the red-blue and the red-green transition zone. 

The gear tooth is subjected to cantilever loading, hence the lower 
end, beyond the loading point remains unstressed i.e. the fringe at the 
outer boundary of the gear tooth is a zeroth fringe order. As we move 
towards the base of the tooth, a source exists at the loading point. 
The stress value increases hence the fringes are ordered in a 
successively increasing order. As we move further down the loading 
point, the stress value decreases. The crack-tip acts as a source and 
in its neighborhood the fringes have to be ordered appropriately. 

The fringes around the crack tip present some unique features. 
When observed in a monochromatic light source, a fringe is seen just 
ahead of the crack tip, along the crack direction <Fig. 5.3). The time 
derivative photographs are taken at loads of 10 Kg, 12 Kg, 15 Kg, 18 
Kg, 20 Kg and 30 Kg respectively. The fringe ahead of the crack 
initially has a fringe order zero which is corroborated by the colour 
code (Fig. 5. 3), As the load is increased, the fringe becomes shorter in 
length. In the photograph obtained using a monochromatic light source, 
the end of the fringe becomes progressively brighter and finally merges 
with the background (Fig. 5-3). In the white light image, it can be seen 
that the fringe has red-blue transition as in fringe order one and 
gradually changes its colour. 

At the crack-tip due to the high stress gradient it is difficult 
to see whether the fringe has red— blue transition or it is black. 
Hence, there is ambiguity in ordering this particular fringe This 
ambiguity is appropriately taken care of in the analysis. Except, this 
fringe, all other fringes are uniquely ordered. 
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Fig. 5.3 Isochromatic fringe patterns of the spur gear tooth with a 
crack at its tensile fillet, for increasing load, under 
monochromatic and white light source. 
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The fringe ordering at the crack-tip is difficult because of two 
reasons. For cracks in structures of finite geometry, a "sink" is 
observed ahead of the crack (Fig. 3. 5). In the present example, apart 
■from this, while going from tensile to compressive fillet the fringe 
order have to go through zero. These two features interact in a complex 
way as the loading is gradually increased. 

5.5 S. I.F evaluation for the crack 

The first step in any photoelastic analysis is calibration of the 
material fringe constant. Using the completely automated procedure 
developed for calibration of a photoelastic material using a circular 
disk under diametral compression, as described in section 3.7 the 
analysis was performed and the material fringe constant was found to be 
14.735 N/mm/fringe- 

Ref erring figure 5.3 it can be seen that the isochromatic loops of 
same order above and below the crack have different sizes, which is a 
distinguishing feature that the loading at the crack is of mixed mode 
in nature. Several factors like the tooth geometry, the position of the 
contact point relative to the crack also affect the isochromatic fringe 
pattern around the crack. The S.I.F analysis is carried out for the 
gear problem for load value of 30 kg using the Atluri and KobayashiC23 
stress field equations as discussed in the previous chapter. 

Seventy five data points were selected, equally distributed to 
represent every geometric feature of the fringe pattern. 

5.5.1 Data processing 

The raw data provided by the data collection software has to be 
processed to transform it from the screen coordinates to real world 
coordinates. This is achieved using the following procedure 
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The aspect ratio oT the MVP-AT screen has to be determined 
initially since the horizontal and vertical pixel scaling is not the 
same- Using the calibration images of a vernier caliper, the 
magnification in terms of pixels/mm is determined for both the X and 
the Y direction. It should be noted that this magnification is not the 
scaling factor of the real world model but scaling factor with respect 
to the pixel size of the screen. For the current settings of the 
monitor the aspect ratio is 0.6818. 

For the purpose of analysis, the positional coordinates of various 
points have to be determined with the crack axis as X axis and the 
^■ack-tip as the origin. 

1. Translation : Depending on the location of the crack tip on the 
screen, the origin has to be shifted from the current screen origin 
( left hand top corner ) to the crack tip. In order to change the 
orientation of the Y axis, the Y coordinates are multiplied by <-l). 



2. Scaling : After the origin is shifted to the crack— tip position, 
the Y coordinates have to be multiplied by the aspect ratio in order 
to obtain the real world distances. 

y„ = 0.6818 Y 

jd. A 

Also the image is a magnified one depending on the type of lens used. 
The crack tip image is recorder using the macro lens. This lens has a 
high magnificat ion and 30 pixels in the X direction represent 1 mm. 


The X and Y coordinates have to be scaled to find the real worl 

coordinates. 

3. Rotation : The crack tip is not located along the positive 


axis. 
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The stress -field equations assume the crack tip to be along the X 
dir ec lion. Hence the data points have to be rotated to align the crack 
tip with the X direction. The angle oT rotation is determined by the 
d-a collection software. In this case the angle was found to be 
28.777^, 

The transformed data available from these transf ormations is 
written to a data file. This file contains the real world X and Y 
coordinates and the fringe order N corresponding to each of the 
selected points. This data is used by the least square software to 
determine the S.I.F parameters- 

5-6 Experl sfierital Results 

The accuracy of any experimental analysis depends upon the 
accuracy of the measurement of the variables. In this case these 
variables are the exact positional coordinates of the data point, the 
location of the origin i.e. the crack-tip, the value of the angle of 
orientation of the crack and the magnification. The accuracy in the 
positional coordinates of the data points and the origin is taken care 
of by the fringe thinning algorithm. A software has been developed to 
measure the angle accurately. It uses human judgment for the purpose of 
Ifjcating the crack direction- Two points, the base and the tip of the 
crack are selected interactively, and the line joining the 2 points and 
extending beyond, is displayed and the 6 value is calculated. 

The least square analysis is then performed on the data collected 
from the processed fringe pattern. The convergence criteria selected is 
the fringe order error minimization- The parameter values are 
resubstituted in the equations and the isochromatic fringe pattern is 
theoretically reconstructed. This is compared with the actual fringe 
pattern by echoing back the data points collected for the analysis. The 
analysis is performed using upto 6 terms of mode I and 6 terms of mode 
II series solution and the results are as shown in Fig, 5. 4. a 
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A sample least square analysis output file is shown in appendix 3. 
Tbe accuracy of the series solution is clearly brought out in the 
appr opr ia te 1 y scaled view of the actual fringe field and the 
reconstructed fringe pattern as shown' in Fig. 5.4.b. The cumulative 
percentage error in the fringe orders v/s the number of terms of the 
series solution is shown in graph in Fig.5.4-c. 

5.7 Closure 

Analysis was carried out for addendum side contact of the gear 
tooth, to find out the optimum number of terms for convergence. The 
criteria for selection of the data points was also discussed. The 
results show that, an higher order solution can faithfully model the 
stress field in this case if appropriate number of terms are used- The 
success of the analysis of the present example provides enough 
confidence to apply photoelasticity to evaluate fracture parameters for 
complex real-life problems. 
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CONCLUSIONS AND SUGGESTIONS FOR FURTHER WORK 


Optical methods in stress analysis have the advantage of 
providing the whole field information of the entire field. A wealth 
of information is provided in the form of fringes. The fringes 
appear in general as bands and the tirst step in analysis is the 
identification of exact fringe location. With the development of 
Digital Image Processing (D.I.P) techniques, sophisticated 
algorithms are available to extract the fringe skeleton. 
Quantitative information can be obtained only if the positional 
coordinates of various points along with their corresponding fringe 
orders are known. The critical problem of the ordering of the 
fringes and the data collection effectively remained a manual one- 
An attempt has been made in this thesis to overcome this bottleneck 
and have complete automation for the photoelastic experimental 
analysis. 

Fringe skeletons are extracted using the algorithm of Ramesh 
and Pramod- The software that was implemented in the PIP environment 
by Ramesh and Pramod was converted to MVP-AT environment in the 
course of this thesis. The fringe thinning algorithm scans the image 
along particular directions and marks all the fringe points. However 
it cannot recognize a particular fringe as a continuous one since, 
in the computer represen tation , the fringes are just a group of 
points without any particular ordering. The first step towards 
automatic fringe ordering lies in the fringe recognition algorithm. 

This is achieved by developing appropriate fringe tracing 
algorithm. Three fringe tracing algorithms are developed, they are 
Simple fringe tracing algorithm, Fringe tracing algorithm with 
directional scan and the complex fringe tracing algorithm. The 
fringe tracing algorithm performs the task of identifying pixels 
forming a particular fringe. The basic methodology is, if a frxnge 
point is selected, the algorithm must be in a position to identify 
^11 contiguous points forming a fringe. The complex fringe tracing 
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elgorithm has been found to work the best for tracing out the fringe 
patterns for test images as well as fringes recorded by actual 

ex per i men ts • 

Fringe tracing algorithm helps to identify and order complete 
fringes if the fringe order of one of the points in the fringe is 
known. To completely order the fringe field, one needs needs to know 
the fringe order for one point each, on different fringes. This was 
made possible by identifying the fringe gradient direction. 

The use of intensity information to identify the fringe 
gradient is established in this thesis. The two dimensional and the 
three dimensional plots revealed that the intensity variation is not 
trivial. Further work is needed to fully utilise this data for 
fringe ordering. 


In the present work collection of data from a the circular disk 
under diametral compression has been automated. This assumes 
significance since, most laboratories use the circular disk under 
diametral compression as a model to evaluate the material fringe 
constant. A semi automatic procedure has been implemented for the 
mixed mode S.I.F analysis for a crack problem. These were made 
possible since the nature of the fringe field is known apriori and 
suitable scan paths were selected for fringe ordering. However, for 
a general class of problems, identification of suitable scan paths 
will be a difficult task- It should be stressed that the problem of 
automation of fringe ordering is very difficult to fit in a standard 
"if-then-else*‘ type of computer algorithm, and one needs 
"intelligent program" to solve it. The future work can be in the 
form of developing an expert system for fringe recognition and 
ordering . 


Fringe ordering is only an intermediate step and one needs 
methodologies to extract field parameters using the fringe data. The 
methodology to extract the stress field parameters in the 
neighborhood of a crack tip from the information of the fringe order 
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and positional coordinates of various ooints 4 - ^ 

pointE was implemented and 

tested for tl.e mode I as well as the mixed mode problems. 


The use of equations proposed by Atluri and KobayashiC23 
for evaluating these stress fieio -=.rameters has been presented. 
Various typographical errors in the reported equations of Atluri and 
KobayashiC23 is corrected in this thesis. 


Gear is one of the most important and most widely used machine 
element. The gear tooth failure commonly occurs due to fracture 
initiated due to a fatigue crack at the tooth fillet. The near tooth 
with a crack has been analysed using photcelasticity . The crack tip 
stress field in the case of a spur gear presents one of the most 
complex fringe patterns for the analysis. Apart from the finite 
geometry of the structure the crack is located in a complex stress 
gradient zone. A semi automatic methodology for ordering the -fringes 
around the crack tip. 

The use of Atluri and Kobayashi's series solution for the 
evaluation of the fracture parameters was discussed- Convergence of 
the series solution depends on a number of factors including the 
choice of the data points, the initial estimates for various 
parameters and the zone of data collection. Analysis was carried out 
for addendum side contact of the gear tooth, to find out the optimum 
number of terms for convergence. The criteria for selection of the 
data points was also discussed. The results show that, an higher 
order solution can faithfully model the stress field in this case if 
appropriate number of terms are used. The success of the analysis of 
the present example provides enough confidence to apply 
photoeiasticity to evaluate fracture parameters for complex 
real-life problems. Further work in the form of analysis of the 
fringe patterns for various gear position to investigate the change 
in the mode I and mode II parameters over the complete contact 
length can be taken up for complete analysis. 
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APPENDIX-1 

MODIFICATIONS REQUIRED TO CONVERT SOFTWARE FROM PIP TO MVP-AT 


One of the major tasks completed in the course of the Dre=^ent 
anvestigation is the conversion of various software modules 
developed around Professional Image Processing ( PIP ) system over 
the last four years to the newly acquired MVP-AT system. Though, 
both PIP and MVP-AT were supplied by the same manufacturer, they 
have not provided appropriate software to convert the programs 
developed from one system to the another. PIP is basically a 
monochrome system on the other hand MVP-AT can support colour image 
processing. In view of the enhanced capabilities of the MVP-AT, even 
the basic structure of the software has to be modified. 

In the following paragraphs, the basic system configuration of 
PIP and MVP-AT are discussed initially. The command structure 
necessary to replace the commands of PIP are summarised in the 
following categories- They are initialisation commands, Input Output 
commands, Point to Point commands, Graphics commands, Statistical 
functions and look up table operations. Finally for various software 
modules, the sequence of quadrants used by it in MVP-AT for 
intermediate operations is summarised. This information is essential 
if one writes a software involving these modules to perform various 


A. 1. System Conf 1 gui^ation 
(i) PIP 

It is a plugin black and white video frame grabber digitizer board 
for the IBM-PC which supports the PIP software- This software has 
libraries for frame buffer I/O, Look Up Table manipulations, 
graphics and imaging functions. The system has a 1Mb frame buffer on 
a PC-8088 machine. The computer languages supported by the PIP 
library are Microsoft— C ( version 4.00 ) and Microsoft FORTRAN 
(version 3.31). 


Description of system parameters : 

1- Frame Buffer : 

An area of PIP memory where related pixel data, which comprises of 
one pixel frame, is stored. The video RAM of the PIP can be 
configured in the following manner: 

a. 4 X 512 X 512 X 8-bit frame buffers. This configuration is 

for a black and white image with the grey values defined by ^ i s 
for every pixel, (quadrant 0 to 3, x snd y coordinates o u x 

b. 1 X 1024 X 0124 X 8-bit frame buffer. ( B to 1023 ) 
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ill) MVf 

ThP MVf'--Al hardware* can be prDgramme>d through the IMABER-Al eo-ftware 
library. The software contains primitives that provide control over 
the frame buffer It), Loo)' Up Table, Graphics and the Imaging 
functions. The system has a 1Mb frame buffer on a PC-386 machine. 
The computer language directly supported by IMAGER-AT is Microsoft-C 

(Version 4 . 00 ) . 

Ill icription of the system parameters : 

1. Frame Buffer : 

The video RAM of the MVP can be configured in the following manner: 

a. 4 X 512 X 512 X 8-bit frame buffers. This configuration is used 

for a blac);. and white image with the grey values defined by 8 

bits for every pixel. ( FB0 to FB3, 0 to 511 each). 

b. 3 X 512 X 512 X 8-bit frame buffers. This configuration is used 

to store coloured images. RED in FB0, GREEN in FBI, BLUE in FB3. The 

frame buffer FB2 is not used. ( FB6, 0 to 511 >. 


c. 2 X 512 X 512 X 16-bit frame buffers. < FB5 and FB6,0 to 511 each) 

d. 1 X 1024 X 0124 X 8-bit frame buffer. ( FB7, 0 to 1023 ) 


A. 2. Modifications necessary to transfer software from PIP to the 
MVP-AT system. 

The PIP software has absolute coordinates referred i.e. 
(0,0)-( 1023, 1023) . Whereas the MVP-AT system is required to be 
configured as 4 X 512 X 512 X 8— bit quadrants. Hence the coordinates 
of a point in any quadrant correspond to the local coordinates i.e. 
(0,0)-(51 1 ,51 1 ) . as shown in Fig.A.l. 


A. 2. 1 . Initialization commands 


PIP : 


f g_ini f mt < baseadd r , mode , speed , c 1 ass , vid-type , zoom ) ; 
fg_chan ( 2 ) ; 
f g_quadm ( 1 ) ; 

■fQ—Pan (0) ; 
fg_scrol 1(0); 
fg_zoom(0) ; 
fg_iniluts( ) ; 
tg_setwin (0,0,511 ,511 ) ; 


The suffix '’fg_" is required for every PIP library command to 
distinguish it from the C library commands. 


the PIP board to a known state. 


The command 


tg_inifmt‘' initializes 



FB0 

FBI 

FB2 

FB3 


FB 

0 

FB2 

FBI 

FB3 
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baseeddr : This parameter speci-fies the o-f-fset o-f the 10 address 
used. Value tor tlie current contiguration : 620. 

mode : Specifies whether the PIP card is configured to allow for 
room resolution. Value tor the current configuration : 1 (allowed). 

speed : Clock speed determining the type of PIP board installed. 

Value tor the current configuration : 0 ( PIP-1024 installed ). 

class : Specifies the video format. 

Value tor the current system configuration : O (512 X 512 EUROPEAN). 

vidtype : Selection between European or American video format. 

Value for the current system configuration : 1 ( EUROPEAN ). 

zoom ; Specifies the resolution. 

Value for the current system configuration : 0 ( full resolution ). 
The current settings are : fg_inifmt<620, 1 ,0,0, 1 ,0) ; 

The command " f g_c han ( chan ) " is used to select the input channel, 
where "chan" is : 

0 : Channel 0, i.e. top camera IP on the PIP card. 

1 : Channel 1, i.e. middle camera IP on the PIP card. 

2 ; Channel 2, i.e. bottom camera IP on the PIP card. 

Value for the current system configuration : 2 < bottom card ). 

The command " f g_quadm ( p) " is used to select the quadrant mode, 
where "p" is : 

0 : 1024 X 1024 mode. 

1 : 512 X 512 X 4 mode. 

Value for the current system configuration : 1 ( 4 quadrants ) . 

The command "fg_pan(x) and fg_scrDl 1 (y ) " is used to set the absolute 
origin position, where <x,y) is the new screen origin. 

Value for the current system configuration : (3,0 ( Top-Left corner 
is the origin ) - 

The command "fg_zoom(r)" is used to select the zoom resolution, 
where "r*‘ is : 

0 (std resolution) - 

1 (compressed image). 

The command f g_ini lute** loads all look up tables with grey scale 

maps • 

The command "fg_setwin" specifies the lower left and upper right 
corners of the working space in all quadrants- „ „ 

Value for the current system configuration ; 0,0,ull,ull 
quad ) . 

MVP-AT : 


if (init(MEMBASE, lOEASE) != OK) < 
exit ( INIT_ERR) ; 

opmode ( PROC, FB0) ; 



dis-f or mat (FfxEO?®, INTLACE , EUROPEAN) ; 
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Di.i t pa t ( 0 , ~ 3 ,0,0) ; 
sync (INTERN, NONE); 
video < ENABLE , EBLir ) ; 
inmode ( MONO ) ; 


HVP-AT is provided with a tile **im2to3.h'‘ which recognizes the MVP 
library tunctions and prefixes them with "im . 


The command ”init*' initialises the board to a default state. This 
command must be used at the begining of the program before any 
c, '’’s to the MVP library. 


The parameter MEMBASE specifies the 64E.; memory segment of the 
board with respect to the PC memory space. The 64K direct 

access window is located in the lower 1Mb of the PC-386 memory 
space- This MVP-AT reserved memory is fully supported by the 
IMAGER-”AT library and the IBM DOS- This allows all the MVP 
frame buffers to be accessable by the bus 64Kb at a time. 

The parameter lOBASE specifies location of board in the PC-Bus I/O 
space . 

The command *'opmode ( mode , quad ) *' selects the operating mode- Every 
command in the IMAGER-AT has to be called from its respective 
operating mode- 

There are four opmodes, 10, GRAPH, PROC, CGRAB for input 
output, graphics routines, imagimg procedures and for functions in 
the continuous grab modes- In the MVP-AT environment, proper opmodes 
need to be set before any library command is used otherwise the 
system behaves unpredictably and may even hang. 

The PIP system does not have any opmodes and once initial ized , the 
system can perform any operation - 

The command '*disf ormat ( f , i , t ) *‘ customizes the display format- 
where f: frequency, i : interlacing parameter, t : format. 

FREQ20 : Selects the frequency of the pixel clock as lOMHz. (Defined 
in the header files ) 

INTLACE : Selects the interlaced monitor. 

EUROPEAN : Selects the european format. 


In PIP this is set by the “fg__inifmt" command itself. 


The command *'outpath( backbuf ,winbuf ,overlaybuf , source) 

select the output display quadrant. Where : 
backbuf : It selects the background frame buffer, 

winbuf : The window frame buffer. 

overlaybuf : Overlay source frame buffer. < Ignored 

negative ) / t h 

source : Selects the nibble in overlaybuf. gnore 

negative ) 


is used to 


when winbuf 
when winbuf 


This command besides displaying the working quadrant can thus be 
used for generating overlays and other Wj.n ow opera io 
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1 hp conrimand sync ( SDL»r t e ) is used Id select the source of the sync 

fj j Q r > bI m 


The source has to be cfianged to *'e>: terna 1 when grabbing an image 
from the camera and changed back when the operation is over. PIP has 
‘'iij sync (source) command for the same purpose and works in a 
similar manner. 

1C command “ vi deo ( di spl ay , source ) ” is used to select the output 
of the frame buffers or the camera on the display, 
display : 1 Enable, 0 Disable, 

source : 0 Camera, 1 Frame Buffer. 

The command *’ inmode ( mode ) *’ is used to select the input data path- 
The MVP system accepts input from cameras of various video 
standards including the separate RGB input and the NTSC. 

MONO: 0 disable, 1 B/W, 2 RGB, 3 NTSC. The PIP does not have such a 
command as it does not support colour images. 

PIP command for changing quadrant : 

f g_dquad ( quad ) ; This displays the quadrant specified by the ’'quad”. 
In HVP-AT one has to use the outpath command. 

outpath< 05 -l ,0,0) ; This selects the quadrant 0 as the display 
quadrant . 


A. 2. 2. I/O 

1. To load the image file from disk to the frame buffer- 
PIP : 

f g_f rdisk ( bsize 5 quad , f name , wbuf ,SBg ) ; 

bsize: Intermediate buffer size, required by DOS, common value 4096 

quad : quadrant. 

fname: filename string. 

wbuf : The PIP working buffer address, common value 5200. 
seg : Portion of work buffer in use, value — 1 - 


MVP-AT : 

opmode < 10 , quad ) ; 
f romdisk ( fname ) ; 

The opmode has to be set to 10 in the required quadrant, 
fname : filename string- 


2- To save the image file to the disk- 


PIP : 

fg todisk ( bsize, quad , fname, wbuf 5seg> ; act'QA 

-r JL L ri3ouired by DOS, common value 

bsize: Intermediate buffer si^^e, require y 

quad : quadrant. 

fname: filename string. t=:r>{ 7 im 

wbuf : The PIP working buffer address, common value u.. 

seg : Portion of work buffer in use, value 
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nVF'-Al : 

opmodE’ ( J 0 ,quad ) ; 
todisf' < 4 namE ) ; 

The opmodr has to bp set to 10 in the required quadrant, 
-fname : filename string. 


3 . Pixel Read/Write : 


PIP : 

■f g p ixr(X||iy) « 

fg_pixw( X ,y , index ) ; 

X f ys They are the x and y coordinates o-f the pixel to be read fron 
the frame buffer ( 0 to 1023 ). 

index : It is the intensity value, an integer variable from 0 to 255 
nyp-AT : 

opmodeC 10, quad ) ; 
pixr ( X ,y ) ; 
pixw ( X , y , pi K ) ; 

>: , y: They are the x and y coordinates of the pixel to be read from 
the frame buffer ( 0 to 511 ). 

pix : An unsigned long variable which accepts the value as red byte 
green byte and last (MSB) blue byte for the colour mode and as the 
lowest byte/bytes for the other modes. 


4. Window to/from buffer Carrayl command. 

PIP s 

f g_winr ( wbuf , seg ) ; 
fg_winw( wbuf ,seg ) ; 

These commands are used to copy the contents of the current display 
window to a buffer or from a buffer. These commands are used for 
storing the window area when it is overlapped by the pop”-up window 
and copying back from the buffer to erase the window. 


MVP-AT : 

These commands are not available in the Mv'P—AT , however these have 
been defined in the header file "curutil 1 . h*^ as follows. 

void Winr<int xl,int yl,int x2,int y2) 
i int i 1 , 1 5 pi ; 
opmode ( 10, quad ) ; 
t=0; 

for ( i=y 1 ; i <= y2; ) f 

rowr ( X 1 5 i 1 5 abs ( x2-x 1 ) , a temp ) ; 
for< pl=0; pi <= abs<x2~xl);++pl)f 
abuf [t3Epl3=atempCpl3; 

> 

t +=1; > 

T 

J' 


void WinwCint xi,int yl^int K2,int y2> 

/ 

t 
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3 n t 1 1 . t , pi; 

opmotlt’ ( 1 0 , quad > ; 
t=0; 


■for(al= yl;3l'> y2’;al+^){ 

for(pl=0; pi <- abs<x2->:l ) ;++pl ) { 
atempC pi D-abu-f [ 1 3 [ pi 3 ; } 
t +=1; 

r oww ( X 1 , i 1 , ( x2-x 1 ) , a temp ) ; 

Ti 

J 

n 

j 

where "abuf" is an array o-f long unsigned 250 X 250 elements, the 
largest array size supported by the compiler and atemp is array of 
unsigned long of 250 elements. It can be seen that this will be 
able to operate with a window size of 250 X 250 of the 512 X 512 
quadrant space- 

However the MVP—AT provides a much more efficient command for 
generating the pop-up menu and such other functions. 

1 l"fi s can be done by using the "drawmode*’ command. 

It sets the pixel replacement policy for the subsequent graphics 
commands - 

Dpmode ( GRAPH , quad ) ; 
drawmode ( 2 ) ; 

This then XOR"s the pixel locations to be overwritten with the new 
values, displays the image when the back ground has to be restored, 
it again XOR wirtes the pixel locations, erasing the current 
window. This is a standard practice used for generating fast 
displays and animations. 

PIP also has a quadrant to quadrant window copy command. 

In MvP—AT one has to use the command “pixblt” to achieve an 
equivalence of *'cpws*'. ''cpws*' is used in various software modules 
repeatedly. The identification of pixblt to do this operation took 
considerable time as there is no semblance of comparison between the 
command names. This command can thus be used directly in place of 
*' c pws ( quad 1 , quad2 ) *' the copy workspace command in PIP? quad! is the 
source quadrant and quad2 is the destination quadrant. 

opmode ( PROC , quad ) ; 

pixbl t ( sbuf ,xl ,yl,x2,y2,dbuf, xSjyS, sdir ,ddir ) ; 
sbuf , dbuf : Source and destination quadrants. 

<xl,yl) , <x2,y2) : Source window coordinates. 

(x3,y3) : Target window starting coordinates. ^ 

sdir , ddir : Source and destination scan direction. Values 
Selecting different source and destination scan directions can 
enable quick rotation of the image through 90, 180, ^/0 degrees. 


A. 2:. 3. Polr*t t o polijt c ormYiarKls 
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5- Logical opr^rators for the image : 


PIP : 

It dops not have any direct logical operators and the corresponding 
pixels have to be read from the two quadrants and then operated upon 
using the functions in Microsof t-C. 

e.g. If two windows have to be logically ANDed then following 
function is required : 

for (i=win>:yC03;i<winKyC2Ii + l;++i) f 

for ( j=win>;y C 1 3 ; j <win>ty Il 33 + 1 j ) { 

c2=f g_pixr <512+i ,512+ j ) ; 
if < <cl != 0) <c2 != 0) ) 

f g_pi>tw( i , j ,250) ; 

> 

\ 

j 

This is replaced by a single command in the MVP-AT system. 

The " in terimage" command performs specified operations between two 
images . 

opmode ( PROC , quad ) ; 

in ter image < sbuf 1 ,sbuf 2 ,dbuf , fn > ; 

sbufl, sbuf2 : The source frame buffers. 

dbuf s The frame buffer in which the resultant image is to be stored 
fn : Specifies the logical operation to be performed. 

0 : ADD 

1 : SUB 

2 : MIN 

3 : MAX 

6 : AND 

7 : XOR 

B : OR . . 

This command is very fast and operates at 1/4 the real time rate. 


A. ^.4. Gr aphics coiinmarids 


6 . Text output on the screen. 


PIP : 

fQ_text( string, size) ; 

fg_tchar (char , size) ; . 4 . - 

string : The text to be output is stored in this string 

size : Specifies the size on the screen, 
char : A single character. 

Writes the string using the current drawing index at 
pen location, e.g. 


the 


current 


void f rn f unc ( ) 



jtit J ; 
t tift r c ; 
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■for ( i =0 ; i <3; 44 i ) 

•f r n buf [ i 3=0; 

■f y_spla nd ( ) ; 

-fQ_rPctf (04qd2.x., 49]+qd2.y, 1824qd2.>:, 


51 1 4qd2. y ) ; 


■fg_sptind (0) ; 

-fg_mDveto(24qd2.K, 5104qd2.y) ; 

-fg_te>;t( "FRINGE: " ,2) ; 

j =0 ; 

do i 

c=getch( ) ; 

if ( <(c >='0' M: c <='9') !! (c =='.')) 

i < 3.) C 

f rn_buf C i 3=c ; 

'^9--^c:har <c ,2) ; /* This function directly 

outputs the number as a character t/ 

i-++; 

j 

y while (c!= crlf && i< 3); 
frn_flag =0; 
f f 1 ush ( stdin ) ; 

j 

This function is used to print FRI^4GE on the screen, then read the 

fringe value, display it on the screen and erase it when comming out 

of the loop- 

MVP-AT : 

opmode ( GRAPH , quad ) ; 

Q tex t < si ze,cnt, string) ; 

size : 1 Normal, 2 Large- 

cnt : Number of bytes to be output. 

string : The text string to be output. 

The character size in case of the MVP~AT is different than the 

cliaracter size in PIP. This command also needs the number of 

bytes in the string, this therefore requires predetermined maximum 
size of the character array 
to be output- e.g. 


Vv^id f rn_f unc ( ) 

r 

t. 

int i ,cl ; 

char c , f rn_tempC 103 ; 

for ( i=0 ; i<3; 4-4-i ) f rn_buf Ei3=0; 
opmode ( GRAPH , 0 ) ; 
drawmode (2 ) ; 
move( 10 , 440 ’ ) ; 
gtext( 1,7, "FRINGE " ) ; 

i=0; 
do f 

c=getch ( ) ; 
c l=toascii (c ) ; 

if ( ((c >='0' &•?< c <='9') !1 <c 


) ) && 


3> 


f rn buf C i 3=c 1 ; 
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operations 


/$ Sance tc har -function is not i 

to bf- per -formed : »/ so -following 

gcvt((cl-4e),2,frn_temp>; 

/♦ Ifus -function is used to convert + ,0 + 

array whose max size is 3 bytes */ variable to a string 

gtext < 1 ,3. ■frn_temp) ; 

/* Ttien tf.is string is displayed »/ 

i -f 4 ; 


3 while (cf= crl-f «;«, i< 3 ). 
move ( 10 j, 440 ) ; ' 

gtext ( 1 ,7, "FRINSE " ) ; 
gtext ( 1 ,3, f rn_temp ) ; 
drawmode ( 0 ) ; 
trn_-flag = 0 ; 
f f lushfstdin) ; 


A- 2 - 5- Statistical functlosis and look up table operations 
The contrast stretching procedure : 

PIP : 

from the interpreter : 

average 00 /t SPATIAL AVERAGING »/ 
histo 5200 /» HISTOGRAM */ 

dq 1 

setindex 0 
cl 0 7 

setindex 255 

dhisto 7.1 10 962 400 255 255 5200 /* Draw the histogram */ 

pause 

dq 0 

modhist 0 %1 72 0 5200 5200 /t Modify the histogram to new range 

*/ 

dq 0 

map 0 0 5200 /t Use that histogram function to map the pixel to new 

intensity values */ 

MVP-AT : 

from the interpreter : 

opmode 2 0 
outpath 0-100 

scaling 20 0 200 255 200 

This command recomputes the intensity vaules between 20 to 200 
and maps then as 0 to 255 linearly. Thus we "stretch" the intensity 
range to the zone of interest only. The inormation contained in 
pixels having intensity values upto 20 and beyond 200 is ignored. */ 

cwb 200 200 -1 15 

This command is used to convert the newly computed values in the 
required format. */ 



ice 


/* thec^P two commands are repeated with difterent 
, 1 , , idang upon the individual image »/ 


intensi ty 


i imi ts 


olutlay 0 200 200 200 000 

/* ItiiE- command actually maps the 
on the screen to the new range »/ 


intensity values of all 


the pixels 


A. 3. Inf oi ttiat Ion stored In different quadrants for 
/ software nksdules. 


va ! 1 o us f unc 1 1 oiis 


1. window_curser < ) ; 


/* Used to select the working window space 
Source tile ; curutill.h t/ 


Input : Original image in quad 0, 

Steps : 

1. Copy original image to quad 2. 

2. Select window and show the clipped part in quad 0. 

not satistactory then copy original image 
from quad 2 to quad 0 and repeat. 

Output : 1. The window portion ot the image in quad 0. 

2. Original image quad 2. 


2. global_threshDld( ) ; /» Used tor thresholding the image. 

Source tile : thresh. h $/ 


Input : Original image in quad 0. 

btfcpZ s 

1. Copy original image to quad 1. 

2. Threshold the image in quad 0. 

3» If it is not satisfactory then copy original image 
from quad 1 to quad 0 and repeat. 


Output 


1 . The thresholded image in quad 0- 
2. Original image quad 1. 


3. hybr id_thin ( ) ; /t For fringe skeleton extraction. 

Source file : orthol.h. */ 

Input ; Thresholded image in quad 0 and Original image in quad 1. 

Steps ; 1 . Row Scan : Scan the image rowwise from quad 0 and write 

the output to quad 2. 

2. Col Scan : Scan the image columnwise troro quad 0 and 

write the output to quad 2. 

3. Di ag_down_thin : Scan the image at Ic-S deg. -from qua^ 0 

and write the output to quad 3. 
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Output 

PROGRAM 
Steps : 


PROGRAM 
Steps : 


uivUun = S.:i,r, thr « dE.g . -from qu.d P- 

and wrjtr the output to quad 3. 

G. AND quadrants 2 and 3 and write the output to quad 2. 

6. OR quadrants 2 and 1 and write the output to quad 0. 

Original image in quad 1 and fringe skeleton superimposed 
on the original image in quadrant 0. 


CALIBl .C 

1 . Load the image in quad 0. 

2. Window_curser ( ) ; 

Input : Original image in quad 0. 

Output : 1, The window portion of the image in quad 0. 

2. Original image quad 2. 

3. 61obal__threshold( ) ; 

Input : Original image in quad 0. 

Output : 1 . The thresholded image in quad 0. 

2. Original image quad 1( window ). 

4. cpws(2,l); /* copy quad 2 to 1 t/ 

5. hybrid_thin ( ) ; 

Input : Thresholded image in quad 0 and Original 
image in quad 1 . 

Output : Original image in quad 1 and fringe 
skeleton superimposed on the original 
image in quadrant 0. 


TILE : 

1. Load the image in quad 0 and save to disk as skn.img. 

2. Load skn.img from disk in quad 0. 

3. Window_curser( ) ; 

Input : Original image in quad 0. 

Output : 1 - The window portion of the image in quad 0 

2. Original image quad 2. 

4. GlDbal_threshold( ) ; 

Input : Original image in quad 0. 

Output : 1. The thresholded image in quad 0. 

2. Original image quad I( window ). 


5. cpws(2,l); /* copy quad 2 to 1 */ 


6. hybrid_thin( 
Input 

Output 


Thresholded image in quad 0 and Original 

image in quad 1 . ... 

Original image in quad 1 and fringe 
skeleton superimposed on the original 
image in quadrant 0. 
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7 . 


li satisfactory then save to disk as skn.imq 
r I * [ ) f * at from s 1 6? p 2 . 


otherwise 


Note on convrrting an image -from the TIF to the MV'P-AT Format : 

AFter converting the image, it 15 necessary to set the window sire 
by the iowin command beFore loading the image in the MVP-AT (Flat) 

format . 



APPENDIX-2 


1 1 1 


ZOOM FUNCTION 


The image obtained -from the CCD camera is 512 
size- ZDoming Tunction is implemented in the 
highlight the zone of importance at the expence of 


X 512 pixels in 
present work to 
the other area. 


1. The image is loaded in quadrant zero and a window is selected to 
mark the portion of the image to be roomed. 

2. The user is asked to input the factor by which the image has to 
be zoomed- This factor has to be selected depending on the window 
fize. If the size of the zoomed window exceeds the quadrant 
dimensions ( 512 X 512 ) then the image gets distorted. 


3. The pixels are then read rowwise and corresponding to each pixel, 
n X n pixels having the same intensity value are drawn in quadrant 
1 . 


4. The coordinate values are then advanced by the factor n. 

This continues till all the rows are scanned. 

Depending on the mask selected the image portion can be zoomed by 
any factor or repeated zooming operations can be performed on the 
image to have a high magnification. The factor of magnification can 
be independent along each axis. However a different factor along x 
and y results in distortion of the picture. 


This function is particularly useful in debugging the automatic 
fringe tracing algorithms. Note that the image contains no new 
information, it only shows a visual magnification by the factor. 
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APPENDIX-3 


Least square output file. 


Convergence 

Error : 0.11000000 





i trns 

parameter values 




1 

1 

4.8206B7 

-0.637504 

0.184755 

0 . 051670 



-0.080121 

0.011736 

-0.131808 

0.011350 



-0.02916B 

-0.017154 

0.004268 


2 

1 

4.849699 

-0.621264 

0.152001 

0.050772 



-0.078510 

0.011293 

-0.131517 

0.011323 



-0 . 029386 

-0.016845 

0.003889 


3 

1 

4.843923 

-0.620714 

0.143154 

0.051264 



-0.078678 

0.011212 

-0.133837 

0.013461 



-0.032331 

-0.016353 

0 . 003672 


4 

1 

4.846506 

-0.618379 

0.135602 

0.051343 



-0.078545 

0.011131 

-0.134098 

0.013818 



-0.032861 

-0.016220 

0.003569 


5 

1 

4.846382 

-0.617720 

0.132626 

0.051415 



-0.078525 

0.011099 

-0.134505 

0.014208 



-0 . 033399 

-0.016120 

0.003516 


6 

1 

4.846701 

-0.617242 

0.130846 

0.051447 



-0.078504 

0.011080 

-0.134632 

0.014347 



-0 . 033598 

-0.016078 

0.003489 



Average Values 


pare 1 3 
parC23 
parC33 


4.842316 

-12.622137 

0.146497 
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par [ 4 3 

0 - 0 L 1318 

par f f»3 = 

- 0.07081 A 

par [63 “ 

0.01 1208 

part 73 = 

- 0.133400 

parT83 = 

0.000000 

par[93 =- 

0.013085 

par r 1 0 3 = 

- 0.031791 

parC113 = 

- 0.016462 

pa r [ 1 2 3 = 

0.003734 


ERROR TABLE 



! y 

i N_data 

i ^4_cal 

i A error 

- 0.58 

0.68 

2.00 

1.90 

5.23 

- 0.64 

0.98 

2.00 

1 .93 

3-69 

- 0.52 

1.35 

2.00 

2.01 

““ 0.26 

- 0.21 

1.56 

2.00 

2.06 

- 2.75 

0.44 

1 .63 

2.00 

1.97 

1.54 

0.76 

1.32 

2.00 

1.97 

1 .71 

0.84 

0.91 

2.00 

2.01 

- 0.67 

0.81 

0.72 

2.00 

2.03 

-1 .69 

0.74 

0.52 

2.00 

2.03 

- 1 .65 

0.68 

0.41 

2.00 

2.01 

- 0.58 

0.51 

0.24 

2.00 

2.04 

- 2.21 

0.36 

0.32 

3.00 

3.02 

- 0 . 72 

0.40 

0.53 

3-00 

3.03 

- 0.85 

0.35 

0.71 

3 . 00 

3.03 

- 0.95 

0.27 

0.84 

3.00 

2.99 

0.21 

0.04 

0.91 

3.00 

2.98 

0.59 

- 0.12 

0.82 

3.00 

3.01 

- 0.34 

- 0 . 23 

0.62 

3 . 00 

3.02 

- 0.67 

- 0.23 

0.54 

3.00 

3.06 

- 1.96 

0.03 

0.47 

4.00 

4.05 

- 1.22 

0.20 

0.43 

4.00 

3.89 

2.86 
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- 0 . 

. 7 t * 

- 0.60 

1 .00 

1.18 

“18 

- ^ . 

.04 

- 0 . 84 

1 .00 

1 .23 

“23 


- 0.93 

- 0.99 

1 .00 

1 .14 

- 14.44 

- 1.10 

-1 .24 

1 .00 

0.90 

9.55 

- 0.53 

- 0.69 

2.00 

1 .92 

4.12 

- 0.55 

-1 .07 

2.00 

1 .94 

3. 15 

- 0.48 

- 1.37 

2.00 

1 .97 

1,31 

- 0.40 

- 1.70 

2.00 

1 .98 

1.00 

- 0.28 

-1 .97 

2.00 

2.05 

- 2.72 

- 0.29 

- 2.25 

2.00 

2.01 

- 0.75 

- 0.47 

- 2.86 

2.00 

1 .97 

1.64 

- 0.83 

- 3.31 

2.00 

2.02 

- 0.89 

1.99 

- 1.92 

2.00 

1 .95 

2.57 

1 .59 

- 1.60 

2.00 

2.03 

- 1.28 

1 .34 

- 1.30 

2.00 

2.06 

- 2.82 

1 . 12 

- 0 . 87 

2.00 

1 .95 

2.49 

0.95 

- 0.70 

2.00 

2.00 

- 0.07 

0.73 

- 0.48 

2.00 

2,05 

- 2.60 

0,64 

- 0.35 

2.00 

1.94 

2.78 

0.46 

- 0.41 

3.00 

2.71 

9.70 

0.51 

- 0.59 

3 . 00 

2.81 

6.37 

8.55 

- 0.87 

3.00 

2.83 

5.53 

0.31 

-1 . 16 

3.00 

2.93 

2.48 

0.03 

- 1 .13 

3.00 

2.87 

4.17 

- 0.12 

- 0.97 

3.00 

2.88 

4.01 

- 0.21 

- 0.87 

3.00 

2.81 

6.21 

- 0.25 

- 0.67 

3.00 

2.89 

3.69 

- 0.02 

- 0.59 

4.00 

3.68 

8.03 

0.20 

- 0.60 

4.00 

3.60 

9.90 

0.26 

- 0.43 

4.00 

3.64 

9.09 

0.66 

- 0.18 

1.00 

1 .33 

- 32.93 

0.91 

- 0,29 

1.00 

1 .31 

- 30.63 

1 .39 

- 0.51 

1 .00 

1.18 

- 18.47 

1.73 

- 0.72 

1 .00 

1.15 

- 14,92 

2.09 

- 0.97 

1.00 

1,12 

- 11.52 



1 15 


1 .02 

- 0.01 

1 .00 

0.89 

10.74 

0.71 

- 0.03 

1 .00 

0.98 

2.07 

0.85 

0.23 

1 .00 

1.31 

- 30.64 

0.74 

0.17 

1 .00 

1.29 

- 28.62 

1.20 

0.38 

1.00 

1 .15 

- 15.42 

1.45 

0.55 

1 .00 

1 .07 

- 7.38 

1 . 66 

0.75 

1 .00 

0.98 

1 .88 

4. 00 

- 1.78 

1 .00 

0.97 

3.15 

4.02 

- 1.27 

1 .00 

0.76 

23.72 

4.02 

- 0.91 

1.00 

0.72 

66 

1.76 

- 0.74 

1 .00 

1 .14 

- 14.29 

2. 12 

- 0.98 

1.00 

1.11 

- 10.84 

2.71 

- 1.30 

1-00 

1 .00 

0.40 

3.06 

-1 .55 

1.00 

1.00 

- 0.19 

3.47 

- 1.72 

1.00 

0.98 

2.15 

2.91 

- 2.53 

2.00 

1.86 

6.78 


The absolute average % error == 6.83 
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